Search

According to the literature review, microbial degradation of polyethylene terephthalate by PETases has been detected effective and eco‐friendly. However, the number of microorganisms capable of such feats is limited with some undesirable bioprospecting results. BTA‐hydrolase has been already reported capable of degrading polyethylene terephthalate. Therefore, mutation by in silico site‐directed mutagenesis means to introduce current isomer of PETase for polyethylene terephthalate degradative capability as a better approach to resolve this issue. This study aimed to use in silico site‐directed mutagenesis to convert a carboxylesterase from Archaeoglobus fulgidus to BTA‐hydrolase from Thermobifida fusca by replacing six amino acids in specific locations. This work was followed by molecular docking analysis with polyethylene terephthalate and polypropylene to compare their interactions. The best‐docked enzyme‐substrate complex was further subjected to molecular dynamics simulation to gauge the binding quality of the BTA‐hydrolase, wild‐type and mutant‐carboxylesterase with only polyethylene terephthalate as a substrate. Results of molecular docking revealed lowest binding energy for the wild‐type carboxylesterase‐polypropylene complex (‐7.5 kcal/mol). The root‐mean‐square deviation value was observed stable for BTA‐hydrolase. Meanwhile, root‐mean‐square fluctuation was assessed with higher fluctuation for the mutated residue Lys178. Consequently, the Rg value for BTA‐hydrolase‐ligand complex (∼1.68 nm) was the lowest compared to the mutant and wild‐type carboxylesterase. The collective data conveyed that mutations imparted a minimal change in the ability of the mutant carboxylesterase to bind to polyethylene terephthalate.

The damage evolution process in 316 austenitic stainless steel at the plastic stage were investigated in this paper by using the in-situ EBSD technique and MD simulation in combination. The damage mechanisms of 316 steel are typified by an intergranular and intragranular mixed mode. The in-situ EBSD results showed that the intergranular damage was the defects formation at the GBs, and the twin-GB intersections and the triple GB are the preferred damage sites. While the intragranular damage was manifested by the mixing of lattice orientation change, intragranular flaws appearance, and twin degradation; The MD simulation results manifested that the initial lattice structure changed after the plastic strain, and the continues increase of the atom displacement resulted in the defect formation. The simulation results also verified that the twin-GB intersection site could be a fragile location, the void preferentially appeared here, and the disorder degree of atoms over here were great. In addition, the defects also generated at the triple GB, which indicated that the triple GB could also be a weak location.

Surfaces with microscale roughness can entail dual-scale hierarchical structures such as the recently reported nano/microstructured surfaces produced in the laboratory (Wang et al. Nature2020, 582, 55−5732494077). However, how the dual-scale hierarchical structured surface affects the apparent wetting/dewetting states of a water droplet, and the transitions between the states are still largely unexplored. Here, we report a systematic large-scale molecular dynamics (MD) simulation study on the wetting/dewetting states of water droplets on various dual-scale nano/near-submicrometer structured surfaces. To this end, we devise slab-water/slab-substrate model systems with a variety of dual-scale surface structures and with different degrees of intrinsic wettability (as measured based on the counterpart smooth surface). The dual-scale hierarchical structure can be described as “nanotexture-on-near-submicrometer-hill”. Depending on three prototypical nanotextures, our MD simulations reveal five possible wetting/dewetting states for a water droplet: (i) Cassie state; (ii) infiltrated upper-valley state; (iii) immersed nanotexture-on-hill state; (iv) infiltrated valley state; and (v) Wenzel state. The transitions between these wetting/dewetting states are strongly dependent on the intrinsic wettability (Ein), the initial location of the water droplet, the height of the nanotextures (H1), and the spacing between nanotextures (W1). Notably, Ein–H1 and Ein–W1 diagrams show that regions of rich wetting/dewetting states can be identified, including regions where between one to five states can coexist.

This work aimed at evaluating the inhibitory effect of ten natural bioactive compounds (1–10) as potential inhibitors of SARS-CoV-2-3CL main protease (PDB ID: 6LU7) and SARS-CoV main proteases (PDB IDs: 2GTB and 3TNT) by molecular docking analysis. The inhibitory effect of all studied compounds was studied with compared to some proposed antiviral drugs which currently used in COVID-19 treatment such as chloroquine, hydroxychloroquine, azithromycin, remdesivir, baloxvir, lopinavir, and favipiravir. Homology modeling and sequence alignment was computed to evaluate the similarity between the SARS-CoV-2-3CL main protease and other SARS-CoV receptors. ADMET properties of all studied compounds were computed and reported. Also, molecular dynamic (MD) simulation was performed on the compound which has the highest binding affinity inside 6LU7 obtained from molecular docking analysis to study it is stability inside receptor in explicit water solvent. Based on molecular docking analysis, we found that caulerpin has the highest binding affinity inside all studied receptors compared to other bioactive compounds and studied drugs. Our homology modeling and sequence alignment showed that SARS-CoV main protease (PDB ID: 3TNT) shares high similarity with 3CLpro (96.00%). Also, ADMET properties confirmed that caulerpin obeys Lipinski’s rule and passes ADMET property, which make it a promising compound to act as a new safe natural drug against SARS-CoV-2-3CL main protease. Finally, MD simulation confirmed that the complex formed between caulerpin and 3CLpro is stable in water explicit and had no major effect on the flexibility of the protein throughout the simulations and provided a suitable basis for our study. Also, binding free energy between caulerpin and 6LU7 confirmed the efficacy of the caulerpin molecule against SARS-CoV-2 main protease. So, this study suggested that caulerpin could be used as a potential candidate in COVID-19 treatment.

The heating process for CuAlNi nanoclusters with cuboctahedral structure was studied by molecular dynamics simulation with the embedded atom method. A diffusionless martensitic transformation from the cuboctahedral (CO) to the icosahedral (ICO) structure before the melting point was found in some CuAlNi nanoclusters. The effects of element composition and cluster size on the thermodynamic stability and structural transformation of trimetallic nanoclusters were investigated. The results show that the CO–ICO transformation temperature was size-dependent and occurred only for small clusters. It was also observed that the CO–ICO transformation temperature firstly dropped and then raised as Al concentration increased, which may be strongly related to the release of excess energy. Furthermore, the surface segregation phenomenon was also observed that the lower surface energy components, such as Al and Cu atoms, segregated to the surface of CuAlNi nanoclusters. The changes in the melting points of CuAlNi nanoclusters with different alloy compositions and cluster sizes were also studied. And furthermore, the physical basis of structural transition of CuAlNi nanoclusters was discussed.

Density Functional Theory (DFT) and Quantitative Structure-Activity Relationship (QSAR) studies were performed on four benzimidazoles (compounds 1–4) and two benzothiazoles (compounds 5 and 6), previously synthesized by our group. The compounds were also investigated for their binding affinity and interactions with the SARS-CoV-2 Mpro (PDB ID: 6LU7) and the human angiotensin-converting enzyme 2 (ACE2) receptor (PDB ID: 6 M18) using a molecular docking approach. Compounds 1, 2, and 3 were found to bind with equal affinity to both targets. Compound 1 showed the highest predictive docking scores, and was further subjected to molecular dynamics (MD) simulation to explain protein stability, ligand properties, and protein–ligand interactions. All compounds were assessed for their structural, physico-chemical, pharmacokinetic, and toxicological properties. Our results suggest that the investigated compounds are potential new drug leads to target SARS-CoV-2.

High-throughput virtual screening methods are known to play key roles in drug discovery. However, the computational docking approaches are often blamed for delivering false-positives in the predicted molecules. Here we demonstrate that incorporation of receptor/ligand dynamics in the computational drug discovery approaches can be instrumental in removing false-positives. In particular, we consider a protein with no well-known inhibitor, Peptidyl-tRNA hydrolase (PTH) and show that a combined protocol of ensemble docking and Molecular Dynamics simulation can efficiently screen out a large set of approximately 1600 candidate molecules. Finally, the addition of metadynamics enhanced sampling technique can assist in probing the resilience of ligands against unbinding and in turn provide a re-ranked selection which can be used for synthetic priority. Overall, the investigation recommends the usage of dynamics in a computational drug discovery program for judging their efficacy.

Overexpression and frequent mutations in FMS-like tyrosine kinase-3 (FLT3) are considered risk factors for severe acute myeloid leukemia (AML). Hyperactive FLT3 induces premature activation of multiple intracellular signaling pathways, resulting in cell proliferation and anti-apoptosis. We conducted the computational modeling studies of 40 pyrimidine-4,6-diamine-based compounds by integrating docking, molecular dynamics, and three-dimensional structure–activity relationship (3D-QSAR). Molecular docking showed that K644, C694, F691, E692, N701, D829, and F830 are critical residues for the binding of ligands at the hydrophobic active site. Molecular dynamics (MD), together with Molecular Mechanics Poison–Boltzmann/Generalized Born Surface Area, i.e., MM-PB(GB)SA, and linear interaction energy (LIE) estimation, provided critical information on the stability and binding affinity of the selected docked compounds. The MD study suggested that the mutation in the gatekeeper residue F691 exhibited a lower binding affinity to the ligand. Although, the mutation in D835 in the activation loop did not exhibit any significant change in the binding energy to the most active compound. We developed the ligand-based comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) models. CoMFA (q2 = 0.802, r2 = 0.983, and QF32 = 0.698) and CoMSIA (q2 = 0.725, r2 = 0.965 and QF32 = 0.668) established the structure–activity relationship (SAR) and showed a reasonable external predictive power. The contour maps from the CoMFA and CoMSIA models could explain valuable information about the favorable and unfavorable positions for chemical group substitution, which can increase or decrease the inhibitory activity of the compounds. In addition, we designed 30 novel compounds, and their predicted pIC50 values were assessed with the CoMSIA model, followed by the assessment of their physicochemical properties, bioavailability, and free energy calculation. The overall outcome could provide valuable information for designing and synthesizing more potent FLT3 inhibitors.

Cardiovascular diseases (CDs) are a major concern in the human race and one of the leading causes of death worldwide. β-Adrenergic receptors (β1-AR and β2-AR) play a crucial role in the overall regulation of cardiac function. In the present study, structure-based virtual screening, machine learning (ML), and a ligand-based similarity search were conducted for the PubChem database against both β1- and β2-AR. Initially, all docked molecules were screened using the threshold binding energy value. Molecules with a better binding affinity were further used for segregation as active and inactive through ML. The pharmacokinetic assessment was carried out on molecules retained in the above step. Further, similarity searching of the ChEMBL and DrugBank databases was performed. From detailed analysis of the above data, four compounds for each of β1- and β2-AR were found to be promising in nature. A number of critical ligand-binding amino acids formed potential hydrogen bonds and hydrophobic interactions. Finally, a molecular dynamics (MD) simulation study of each molecule bound with the respective target was performed. A number of parameters obtained from the MD simulation trajectories were calculated and substantiated the stability between the protein-ligand complex. Hence, it can be postulated that the final molecules might be crucial for CDs subjected to experimental validation.

International audience; Poly(vinylidene fluoride) (PVDF) is identified as a piezoelectric polymer and has attracted a great deal of attention for various applications. It has been proved experimentally that adding Fe3O4 nanoparticles improves β phase contribution, which has the highest piezoelectric coefficient among all the PVDF crystalline forms. Here, we used molecular dynamics simulation to investigate the physical properties and conformational changes of α-PVDF in the presence of magnetite slabs with different surface chemistry. The results show that H-bond interactions between OH groups on the surface of magnetite and PVDF chains stabilize the system and cause the α → β phase transition of the polymer. Moreover, replacing 20% of hydroxyl groups by carboxylic acids reduces the extent of the phase transition. In turn, hexanoic acid grafting induces the highest boosting effect, among all the carboxylic acid modifications of the magnetite, and the interaction of PVDF chains with the magnetite slab increases as the hexanoic acid:OH ratio decreases.

The recognition of carbohydrates by lectins plays key roles in diverse cellular processes such as cellular adhesion, proliferation, and apoptosis, which makes it a therapeutic target of significance against cancers. One of the most functionally active lectins, galectin-3 is distinctively known for its specific binding affinity toward β-galactoside. However, despite the prevalence of high-resolution crystallographic structures, the mechanistic basis and more significantly, the dynamic process underlying carbohydrate recognition by galectin-3 are currently elusive. To this end, we employed extensive Molecular Dynamics simulations to unravel the complete binding event of human galectin-3 with its native natural ligand N-acetyllactosamine (LacNAc) at atomic precision. The simulation trajectory demonstrates that the oligosaccharide diffuses around the protein and eventually identifies and binds to the biologically designated binding site of galectin-3 in real time. The simulated bound pose correlates with the crystallographic pose with atomic-level accuracy and recapitulates the signature stabilizing galectin-3/oligosaccharide interactions. The recognition pathway also reveals a set of transient non-native ligand poses in its course to the receptor. Interestingly, kinetic analysis in combination with a residue-level picture revealed that the key to the efficacy of a more active structural variant of the LacNAc lay in the ligand’s resilience against disassociation from galectin-3. By catching the ligand in the act of finding its target, our investigations elucidate the detailed recognition mechanism of the carbohydrate-binding domain of galectin-3 and underscore the importance of ligand–target binary complex residence time in understanding the structure–activity relationship of cognate ligands.

Amyloid peptides are known to self-assemble into larger aggregates that are linked to the pathogenesis of many neurodegenerative disorders. In contrast to primary nucleation, recent experimental and theoretical studies have shown that many toxic oligomeric species are generated through secondary processes on a pre-existing fibrillar surface. Nucleation, for example, can also occur along the surface of a pre-existing fibril—secondary nucleation—as opposed to the primary one. However, explicit pathways are still not clear. In this study, we use molecular dynamics simulation to explore the free energy landscape of a free Abeta monomer binding to an existing fibrillar surface. We specifically look into several potential Abeta structural precursors that might precede some secondary events, including elongation and secondary nucleation. We find that the overall process of surface-dependent events can be described at least by the following three stages: 1. Free diffusion 2. Downhill guiding 3. Dock and lock. And we show that the outcome of adding a new monomer onto a pre-existing fibril is pathway-dependent, which leads to different secondary processes. To understand structural details, we have identified several monomeric amyloid precursors over the fibrillar surfaces and characterize their heterogeneity using a probability contact map analysis. Using the frustration analysis (a bioinformatics tool), we show that surface heterogeneity correlates with the energy frustration of specific local residues that form binding sites on the fibrillar structure. We further investigate the helical twisting of protofilaments of different sizes and observe a length dependence on the filament twisting. This work presents a comprehensive survey over the properties of fibril growth using a combination of several openMM-based platforms, including the GPU-enabled openAWSEM package for coarse-grained modeling, MDTraj for trajectory analysis, and pyEMMA for free energy calculation. This combined approach makes long-timescale simulation for aggregation systems as well as all-in-one analysis feasible. We show that this protocol allows us to explore fibril stability, surface binding affinity/heterogeneity, as well as fibrillar twisting. All these properties are important for understanding the molecular mechanism of surface-catalyzed secondary processes of fibril growth.

SMYD3 is a SET-domain-containing methyltransferase that catalyzes the transfer of methyl groups onto lysine residues of substrate proteins. Methylation of MAP3K2 by SMYD3 has been implicated in Ras-driven tumorigenesis, which makes SMYD3 a potential target for cancer therapy. Of all SMYD family proteins, SMYD3 adopt a closed conformation in a crystal structure. Several studies have suggested that the conformational changes between the open and closed forms may regulate the catalytic activity of SMYD3. In this work, we carried out extensive molecular dynamics simulations on a series of complexes with a total of 21 μs sampling to investigate the conformational changes of SMYD3 and unveil the molecular mechanisms. Based on the C-terminal domain movements, the simulated models could be depicted in three different conformational states: the closed, intermediate and open states. Only in the case that both the methyl donor binding pocket and the target lysine-binding channel had bound species did the simulations show SMYD3 maintaining its conformation in the closed state, indicative of a synergetic effect of the cofactors and target lysine on regulating the conformational change of SMYD3. In addition, we performed analyses in terms of structure and energy to shed light on how the two regions might regulate the C-terminal domain movement. This mechanistic study provided insights into the relationship between the conformational change and the methyltransferase activity of SMYD3. The more complete understanding of the conformational dynamics developed here together with further work may lay a foundation for the rational drug design of SMYD3 inhibitors.

Fresh two-dimensional molybdenum disulfide (MoS2) absorbs the hydrocarbon contaminations in the ambient air and makes surface aging. To understand how the surface aging influences the interactions between MoS2 and biomolecules is important in the biomedical applications. Here, employing all-atom molecular dynamics simulations, we investigated the interactions of the fresh and aged MoS2 nanosheets with the lipid membranes of different components. Our results demonstrate that both the fresh and aged MoS2 nanosheets can spontaneously insert into the bilayer membranes. However, the fresh MoS2 nanosheet displays significantly stronger interaction and then has a larger penetration depth than the aged counterpart, regardless of the lipid components. The calculations of potential mean forces through the umbrella sampling further confirm that the insertion of fresh MoS2 into the lipid membranes is more energetically favorable. Moreover, we found that the fresh MoS2 nanosheet can cause a larger damage to the integrity of lipid membranes than the aged one. This work provides insightful understandings of the surface-aging-dependent interactions of the MoS2 nanosheets with biomembranes, which could facilitate the design of novel MoS2-based nanodevices with advanced surface properties.

Ion-ion interactions (salt bridges) between favorable pairs of charged residues are important for the conformational stability of proteins. Molecular dynamic (MD) simulations are useful for elucidating the interactions among charged residues fluctuating in solution. However, the quality of MD results depends strongly on the force fields used. In this study, we compared the strengths of salt bridges among force fields by performing MD simulations using the CutA1 protein (trimer) from the hyperthermophile Pyrococcus horikoshii (PhCutA1), which has an unusually large proportion of charged residues. The force fields Chemistry at HARvard Macromolecular Mechanics (Charmm)27, Assisted Model Building and Energy Refinement (Amber)99sb, Amber14sb, GROningen Molecular Simulation (Gromos)43a1, and Gromos53a6 were used in combination with two different water models, tip3p (for Charmm27, Amber99sb, and Amber14sb) and simple point charge/extended (for Amber99sb, Gromos43a1, and Gromos53a6), yielding a total of six combinations. The RMSDs of all Cα atoms of PhCutA1 were similar among force fields, except for Charmm27, during 400-ns MD simulations at 300 K; however, the radius of gyration (Rg ) was greater for Amber99sb and shorter for Gromos43a1. The average strengths of salt bridges for each positively charged residue did not differ greatly among force fields, but the strengths at specific sites within the structure depended sensitively on the force field used. In the case of the Gromos group, positively charged residues could engage in favorable interactions with many more charged residues than in the other force fields, especially in loop regions; consequently, the apparent strength at each site was lower.

Histone deacetylase (HDAC) 8 has been implicated as a potential therapeutic target in a variety of cancers, neurodegenerative disorders, metabolic dysregulation and autoimmune and inflammatory diseases. Several nonselective HDAC inhibitors have been co-crystallized with HDAC8. Molecular dynamics (MD) studies may yield valuable information on the structural stabilities of the complexes over time as determined by various pharmacophore features of the co-crystallized inhibitors. Here, using 11 unmodified X-ray crystal structures of human HDAC8 (complexes) structure-based pharmacophore models were built and clustered based on distance - a function of the number of common pharmacophore features and the root-mean-squared displacement between the matching features. Based on this information, a total of seven complexes (1T64, 1W22, 3RQD, 3SFF, 3F0R, 5VI6 and 5FCW) were submitted to unrestrained 50 ns-MD simulations using nanoscale MD (NAMD) software. 1T64 (HDAC8 in complex with TSA) was found to show the highest stability over time, presumably because of the TSA's ability to span HDAC8 catalytic channel and form a strong ionic interaction with zinc metal ion. Other stable complexes were 1W22, 3SFF, 3F0R and 5FCW. However, 3RQD and 5VI6 showed relative instability over 50 ns time period. This may be attributed to bulkiness of the capping groups of both largazole thiol and trapoxin A, making them unable to fit well into the active site of HDAC8. They rather formed steric clashes with residues on loop regions near the entrance to the channel. Thus, 1T64 and similar crystal structures may be good candidates for HDAC8 structural dynamics studies and inhibitor design.Communicated by Ramaswamy H. Sarma.

Significance For living cells to function, proteins must efficiently navigate the densely packed cytosol. Protein diffusion is slowed down by high viscosity and can come to a complete halt because of nonspecific binding and aggregation. Using molecular dynamics simulations, we develop a detailed description of protein diffusion in concentrated protein solution. We confirm that soluble proteins in concentrated solutions diffuse not as isolated particles, but as members of transient clusters between which they constantly exchange. Nonspecific protein binding and the formation of dynamic clusters nearly quantitatively account for the high viscosity and slow diffusivity in concentrated protein solutions, consistent with the Stokes–Einstein relations., We develop a detailed description of protein translational and rotational diffusion in concentrated solution on the basis of all-atom molecular dynamics simulations in explicit solvent. Our systems contain up to 540 fully flexible proteins with 3.6 million atoms. In concentrated protein solutions (100 mg/mL and higher), the proteins ubiquitin and lysozyme, as well as the protein domains third IgG-binding domain of protein G and villin headpiece, diffuse not as isolated particles, but as members of transient clusters between which they constantly exchange. A dynamic cluster model nearly quantitatively explains the increase in viscosity and the decrease in protein diffusivity with protein volume fraction, which both exceed the predictions from widely used colloid models. The Stokes–Einstein relations for translational and rotational diffusion remain valid, but the effective hydrodynamic radius grows linearly with protein volume fraction. This increase follows the observed increase in cluster size and explains the more dramatic slowdown of protein rotation compared with translation. Baxter’s sticky-sphere model of colloidal suspensions captures the concentration dependence of cluster size, viscosity, and rotational and translational diffusion. The consistency between simulations and experiments for a diverse set of soluble globular proteins indicates that the cluster model applies broadly to concentrated protein solutions, with equilibrium dissociation constants for nonspecific protein–protein binding in the Kd ≈ 10-mM regime.

Joining different parts is one of the crucial components of designing/engineering of materials. Presently, the current energy efficient low weight automotive and aerospace components consist of a different class of materials, such as metals, polymers, ceramics, etc. Joining these components remains a challenge. Here, we demonstrate metal (aluminum) and polymer (Polypropylene, pp) joining using mechanical friction. The detailed characterization clearly demonstrates that atomically locked interfaces are formed in such joining and no chemical bonds are formed during the joining. Also, a waterproof and strong interface is formed in such a process. Fully atomistic molecular dynamics simulations were also carried out in order to further gain insights on the joining process., by Arpan Rout, Prafull Pandey, Eliezer Fernando Oliveira, Pedro Alvesda Silva Autreto, Anurag Gumaste, Amit Singh, Douglas Soares Galvao, Amit Arora and Chandra Sekhar Tiwary

In this study, we provide a time-dependent (td-) mechanical model, taking advantage of molecular dynamics (MD) simulations, quasiharmonic analysis of MD trajectories and td-linear response theories (td-LRT) to describe vibrational energy redistribution within the protein matrix. The theoretical description explains the observed biphasic responses of specific residues in myoglobin to CO-photolysis and photoexcitation on heme. The fast responses are found triggered by impulsive forces and propagated mainly by principal modes -1. The predicted fast responses for individual atoms are then used to study signal propagation within protein matrix and signals are found to propagate ∼ 8 times faster across helices (4076 m/s) than within the helices, suggesting the importance of tertiary packing in proteins’ sensitivity to external perturbations. We further develop a method to integrate multiple intramolecular signal pathways and discover frequent “communicators”. These communicators are found evolutionarily conserved including those distant from the heme.

The combination of molecular dynamics (MD) simulations and small-angle X-ray scattering (SAXS), called the MD-SAXS method, is efficient for investigating protein dynamics. To overcome the time-scale limitation of all-atom MD simulations, coarse-grained (CG) representations are often utilized for biomolecular simulations. In this study, we propose a method to combine CG MD simulations with SAXS, termed the CG-MD-SAXS method. In the CG-MD-SAXS method, the scattering factors of CG particles for proteins and nucleic acids are evaluated using high-resolution structural data in the Protein Data Bank, and the excluded volume and the hydration shell are modeled using two adjustable parameters to incorporate solvent effects. To avoid overfitting, only the two parameters are adjusted for an entire structure ensemble. To verify the developed method, theoretical SAXS profiles for various proteins, DNA/RNA, and a protein-RNA complex are compared with both experimental profiles and theoretical profiles obtained by the all-atom representation. In the present study, we applied the CG-MD-SAXS method to the Swi5-Sfr1 complex and three types of nucleosomes to obtain reliable ensemble models consistent with the experimental SAXS data.

Background: The SARS-CoV-2 main protease (Mpro) is an attractive target in the COVID-19 drug development process. It catalyzes the polyprotein's translation from viral RNA and specifies a particular cleavage site. Due to the absence of identical cleavage specificity in human cell proteases, targeting Mpro with chemical compounds can obstruct the replication of the virus. Methods: To explore the potential binding mechanisms of 1,2,3-triazole scaffolds in comparison to co-crystallized inhibitors 11a and 11b towards Mpro, we herein utilized molecular dynamics and enhanced sampling simulation studies. Results and conclusion: All the 1,2,3-triazole scaffolds interacted with catalytic residues (Cys145 and His41) and binding pocket residues of Mpro involving Met165, Glu166, Ser144, Gln189, His163, and Met49. Furthermore, the adequate binding free energy and potential mean force of the topmost compound 3h was comparable to the experimental inhibitors 11a and 11b of Mpro. Overall, the current analysis could be beneficial in developing the SARS-CoV-2 Mpro potential inhibitors. © 2022 Elsevier Ltd 5058; Department of Biotechnology, Ministry of Science and Technology, India, DBT; Indian Council of Medical Research, ICMR; Council of Scientific and Industrial Research, India, CSIR; Ministry of Education and Science of the Russian Federation, Minobrnauka: 075-15-2020-777 We gratefully acknowledge to the Director, CSIR-Institute of Himalayan Bioresource Technology, Palampur for providing the facilities to carry out this work. This research received no specific grant from any funding agency. The work was carried out under the aegis of the Himalayan Centre for High-throughput Computational Biology (HiCHiCoB), a BIC supported by DBT. The CSIR support in the form of projects MLP:0201 and OLP:0043 for bioinformatics studies is highly acknowledged. R. S. expresses gratitude to the Indian Council of Medical Research, New Delhi, India for providing Senior Research Fellowship. G.V.Z. acknowledge the Ministry of Science and Higher Education of the Russian Federation within the framework of the grant agreement as government subsidies from the Federal budget in accordance with paragraph 4 of article 78.1 of the Budget Code of the Russian Federation (Moscow, October 1, 2020, No. 075-15-2020-777). This manuscript represents CSIR-IHBT communication no. 5058. We gratefully acknowledge to the Director, CSIR-Institute of Himalayan Bioresource Technology, Palampur for providing the facilities to carry out this work. This research received no specific grant from any funding agency. The work was carried out under the aegis of the Himalayan Centre for High-throughput Computational Biology (HiCHiCoB), a BIC supported by DBT. The CSIR support in the form of projects MLP:0201 and OLP:0043 for bioinformatics studies is highly acknowledged. R. S. expresses gratitude to the Indian Council of Medical Research, New Delhi, India for providing Senior Research Fellowship. G.V.Z. acknowledge the Ministry of Science and Higher Education of the Russian Federation within the framework of the grant agreement as government subsidies from the Federal budget in accordance with paragraph 4 of article 78.1 of the Budget Code of the Russian Federation (Moscow, October 1, 2020, No. 075-15-2020-777). This manuscript represents CSIR-IHBT communication no. 5058.

Androgen receptor (AR) is an important therapeutic target for the treatment of diseases such as prostate cancer, hypogonadism, muscle wasting, etc. In this study, the complex structures of the AR ligand-binding domain (LBD) with fifteen ligands were analyzed by molecular dynamics simulations combined with the alanine-scanning-interaction-entropy method (ASIE). The quantitative free energy contributions of the pocket residues were obtained and hotspot residues are quantitatively identified. Our calculation shows that that these hotspot residues are predominantly hydrophobic and their interactions with binding ligands are mainly van der Waals interactions. The total binding free energies obtained by summing over binding contributions by individual residues are in good correlation with the experimental binding data. The current quantitative analysis of binding mechanism of AR to ligands provides important insight on the design of future inhibitors.

In this study, we have investigated the binding mechanism of two FDA approved drugs (ivermectin and levosalbutamol) with the spike protein of SARs-CoV-2 using three different computational modeling techniques. Molecular docking results predict that ivermectin shows a large binding affinity for spike protein (-9.0 kcal/mol) compared to levosalbutamol (-4.1 kcal/mol). Ivermectin binds with GLN492, GLN493, GLY496 and TRY505 residues in the spike protein through hydrogen bonds and levosalbutamol binds with TYR453 and TYR505 residues. Using density functional theory (DFT) studies, we have calculated the binding energies between ivermectin and levosalbutamol with residues in spike protein which favor their binding are − 17.8 kcal/mol and − 20.08 kcal/mol, respectively. The natural bond orbital (NBO) charge analysis has been performed to estimate the amount of charge transfer that occurred by two drugs during interaction with residues. Molecular dynamics (MD) study confirms the stability of spike protein bound with ivermectin through RMSD and RMSF analyses. Three different computer modeling techniques reveal that ivermectin is more stable than levosalbutamol in the active site of spike protein where hACE2 binds. Therefore, ivermectin can be a suitable inhibitor for SARS-CoV-2 to enter into the human cell through hACE2.

Corona Virus Disease 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome coronavirus (SARS CoV-2) has been declared a worldwide pandemic by WHO recently. The complete understanding of the complex genomic structure of SARS CoV-2 has enabled the use of computational tools in search of SARS CoV-2 inhibitors against the multiple proteins responsible for its entry and multiplication in human cells. With this endeavor, 177 natural, anti-viral chemical entities and their derivatives, selected through the critical analysis of the literatures, were studied using pharmacophore screening followed by molecular docking against RNA dependent RNA polymerase and main protease. The identified hits have been subjected to molecular dynamic simulations to study the stability of ligand-protein complexes followed by ADMET analysis and Lipinski filters to confirm their drug likeliness. It has led to an important start point in the drug discovery and development of therapeutic agents against SARS CoV-2., Graphical abstract Image, graphical abstract

Fascaplysin is a natural marine product originating from sponges, attracting widespread attention due to its potential inhibitory activities against CDK4. However, its clinical application has been largely limited because of serious adverse effects caused by planar skeleton. To reduce the serious adverse effects, 18 tetrahydro-β-carboline analogs (compounds 6a-i and 7a-i) were designed and synthesized via breaking the planarity of fascaplysin, and the biological activities of the synthesized compounds were evaluated by MTT assay and CDK4/CycD3 enzyme inhibition assay. The title compounds showed varying degrees of inhibitory activities, especially the cytotoxicity of compound 6c against HeLa cells (IC50 = 1.03 ± 0.19 μM) with quite weak cytotoxicity toward the normal cells WI-38 (IC50 = 311.51 ± 56.06 μM), and the kinase inhibition test indicated that compound 6c was a potential CDK4 inhibitor. In order to further compare the action mechanisms of planar and nonplanar molecules on CDK4, the studied complexes of CDK4 bound with fascaplysin and three representative compounds (compound 6a-c) with bioactivities gradient were constructed by molecular docking and further verified through molecular dynamic simulation, which identified the key residues contributing to the ligands’ binding. By comparing the binding modes of the constructed systems, it could be found that the residues contributing significantly to compound 6c′s binding were highly consistent with those contributing significantly to fascaplysin’s binding. Through the design, synthesis of the nonplanar fascaplysin derivatives, and binding mechanism analysis, some valuable hints for the discovery of antitumor drug candidates could be provided.

Novel coronavirus SARS-CoV-2continues tospread rapidly worldwide and causing serious health and economic loss. In the absence of any effective treatment, various in-silico approaches are being explored towards the therapeutic discovery against COVID-19. Targeting multiple key enzymes of SARS-CoV-2 with a single potential drug could be an important in-silico strategy to tackle the therapeutic emergency. A number of Food and Drug Administration (FDA) approved drugs entered into clinical stages were originated from multi-target approaches with an increased rate, 16-21% between 2015 and 2017. In this study, we selected an FDA-approved library (Prestwick Chemical Library of 1520 compounds) and implemented in-silico virtual screening against multiple protein targets of SARS-CoV-2 on the Glide module of Schrodinger software (release 2020-1). Compounds were analyzed for their docking scores and the top-ranked against each targeted protein were further subjected to Molecular Dynamics (MD) simulations to assess the binding stability of ligand-protein complexes. A multi-targeting approach was optimized that enabled the analysis of several compounds' binding efficiency with more than one protein targets. It was demonstrated that Diosmin (6) showed the highest binding affinity towards multiple targets with binding free energy (kcal/mol) values of -63.39 (nsp3); -62.89 (nsp9); -31.23 (nsp12); and -65.58 (nsp15). Therefore, our results suggests that Diosmin (6) possesses multi-targeting capability, a potent inhibitor of various non-structural proteins of SARS-CoV-2, and thus it deserves further validation experiments before using as a therapeutic against COVID-19 disease.

Recent studies have identified cyclic GMP-AMP synthase (cGAS) as an important target for treating autoimmune diseases, and several inhibitors of human cGAS (hcGAS) and their structures in complexation with hcGAS have been reported. However, the mechanisms via which these inhibitors interact with hcGAS are not completely understood. Here, we aimed to assess the performance of molecular mechanics/Poisson&ndash, Boltzmann solvent-accessible surface area (MM/PBSA) in evaluating the binding affinity of various hcGAS inhibitors and to elucidate their detailed interactions with hcGAS from an energetic viewpoint. Using molecular dynamics (MD) simulation and MM/PBSA approaches, the estimated free energies were in good agreement with the experimental ones, with a Pearson&rsquo, s correlation coefficient and Spearman&rsquo, s rank coefficient of 0.67 and 0.46, respectively. In per-residue energy decomposition analysis, four residues, K362, R376, Y436, and K439 in hcGAS were found to contribute significantly to the binding with inhibitors via hydrogen bonding, salt bridges, and various &pi, interactions, such as &pi, &middot, &pi, stacking, cation&middot, hydroxyl&middot, and alkyl&middot, interactions. In addition, we discussed other key interactions between specific residues and ligands, in particular, between H363 and JUJ, F379 and 9BY, and H437 and 8ZM. The sandwiched structures of the inhibitor bound to the guanidinium group of R376 and the phenyl ring of Y436 were also consistent with the experimental data. The results indicated that MM/PBSA in combination with other virtual screening methods, could be a reliable approach to discover new hcGAS inhibitors and thus is valuable for potential treatments of cGAS-dependent inflammatory diseases.

SARS-CoV-2 has rapidly emerged as a global pandemic with high infection rate. At present, there is no drug available for this deadly disease. Recently, Mpro (Main Protease) enzyme has been identified as essential proteins for the survival of this virus. In the present work, Lipinski's rules and molecular docking have been performed to identify plausible inhibitors of Mpro using food compounds. For virtual screening, a database of food compounds was downloaded and then filtered using Lipinski's rule of five. Then, molecular docking was accomplished to identify hits using Mpro protein as the target enzyme. This led to identification of a Spermidine derivative as a hit. In the next step, Spermidine derivatives were collected from PubMed and screened for their binding with Mpro protein. In addition, molecular dynamic simulations (200 ns) were executed to get additional information. Some of the compounds are found to have strong affinity for Mpro, therefore these hits could be used to develop a therapeutic agent for SARS-CoV-2., Graphical abstract Image 1

In this work we aimed to perform an in silico predictive screening, docking and molecular dynamic study to identify 1,2,3-triazole-phthalimide derivatives as drug candidates against SARS-CoV-2. The in silico prediction of pharmacokinetic and toxicological properties of hundred one 1,2,3-triazole-phtalimide derivatives, obtained from SciFinder® library, were investigated. Compounds that did not show good gastrointestinal absorption, violated the Lipinski's rules, proved to be positive for the AMES test, and showed to be hepatotoxic or immunotoxic in our ADMET analysis, were filtered out of our study. The hit compounds were further subjected to molecular docking on SARS-CoV-2 target proteins. The ADMET analysis revealed that 43 derivatives violated the Lipinski’s rules and 51 other compounds showed to be positive for the toxicity test. Seven 1,2,3-triazole-phthalimide derivatives (A7, A8, B05, E35, E38, E39, and E40) were selected for molecular docking and MFCC—ab initio analysis. The results of molecular docking pointed the derivative E40 as a promising compound interacting with multiple target proteins of SARS-CoV-2. The complex E40-Mpro was found to have minimum binding energy of −10.26 kcal/mol and a general energy balance, calculated by the quantum mechanical analysis, of −8.63 eV. MD simulation and MMGBSA calculations confirmed that the derivatives E38 and E40 have high binding energies of −63.47 ± 3 and −63.31 ± 7 kcal/mol against SARS-CoV-2 main protease. In addition, the derivative E40 exhibited excellent interaction values and inhibitory potential against SAR-Cov-2 main protease and viral nucleocapsid proteins, suggesting this derivative as a potent antiviral for the treatment and/or prophylaxis of COVID-19. Communicated by Ramaswamy H. Sarma

To further explore the long-term aging behaviors of bitumen from the multiscale perspectives, the experimental characterization and molecular dynamics (MD) simulation methods were performed. Series of chemical properties for the virgin and various aged bitumen were evaluated using the TCL-FID, ATR-FTIR, Elemental analysis and GPC tests. The molecular models of virgin and aged bitumen were established firstly, and the influence of long-term aging on the thermodynamics properties was predicted from the MD simulation results. The experimental results revealed that with the aging degree deepened, the resin and asphaltene fractions both increased dramatically, which resulted in the increment of average molecular weight and the more uneven molecular weight distribution in aged bitumen. Moreover, the aging of bitumen led to the increase of the oxidized functional groups (C O and S O) index, oxygen content, aromaticity and polarity, while the carbon, hydrogen element contents and the H/C ratio reduced. The density values from MD simulation agreed well with the experimental results, which significantly validated the reliability of molecular models for the virgin and different aged bitumen binders. The MD simulation results demonstrated that the long-term aging remarkably improved the cohesive energy density, solubility parameter and activation energy, however it deteriorated the surface free energy, work of cohesion and self-diffusion coefficient of the bitumen molecular system. This study develops the molecular models of virgin and aged bitumen with different long-term aging degrees, and provides a fundamental understanding regarding the influence of long-term aging influence on the chemo-physical and thermodynamic properties of bitumen.

International audience; A combined molecular dynamics (MD) simulation and density functional theory (DFT) calculations were used to investigate the effect of the acetonitrile (ACN) solution on the absorption and the charge transfer in a D205 indoline dye at the interface of a TiO2 nanoparticle. DFT calculations were carried out to estimate the equilibrium geometry of a small Ti30O62H4-anatase nanoparticle and to derive interaction parameters for bidentate binding of a D205 dye to TiO2. A series of the D205 dye configurations anchored to Ti30O62H4 were generated from the MD simulations and used as input for the time-dependent DFT (TD-DFT) calculations. We found that the immersion of the D205-Ti30O62H4 composite system into the polar ACN environment results in conformational changes of the dye moiety, which are accompanied by the solvatochromic redshift of its long-wavelength absorption band up to 1900 cm-1. Our results show that the HOMO and LUMO energy level alignment of the dye and the nanoparticle suggests a favourable driving force for charge transfer from the dye to TiO2. The MD/TD-DFT-estimated solvatochromic shift for the D205 dye at the TiO2/ACN interface is in good agreement with the experiment, showing that such computational approach enables reliable predictions of optical properties of other dye-TiO2 composites in solution

Abstract Membrane proteins govern critical cellular processes and are central to human health and associated disease. Understanding of membrane protein function is obscured by the vast ranges of structural dynamics—both in the spatial and time regime—displayed in the protein and surrounding membrane. The membrane lipids have emerged as allosteric modulators of membrane protein function, which further adds to the complexity. In this review, we discuss several examples of membrane dependency. A particular focus is on how molecular dynamics (MD) simulation have aided to map membrane protein dynamics and how enhanced sampling methods can enable observing the otherwise inaccessible biological time scale. Also, time-resolved X-ray scattering in solution is highlighted as a powerful tool to track membrane protein dynamics, in particular when combined with MD simulation to identify transient intermediate states. Finally, we discuss future directions of how to further develop this promising approach to determine structural dynamics of both the protein and the surrounding lipids. Graphic Abstract

Cardiac muscle thin filaments are composed of actin, tropomyosin, and troponin that change conformation in response to Ca2+ binding, triggering muscle contraction. Human cardiac troponin C (cTnC) is the Ca2+-sensing component of the thin filament. It contains structural sites (III/IV) that bind both Ca2+ and Mg2+ and a regulatory site (II) that has been thought to bind only Ca2+. Binding of Ca2+ at this site initiates a series of conformational changes that culminate in force production. However, the mechanisms that underpin the regulation of binding at site II remain unclear. Here, we have quantified the interaction between site II and Ca2+/Mg2+ through isothermal titration calorimetry and thermodynamic integration simulations. Direct and competitive binding titrations with WT N-terminal cTnC and full-length cTnC indicate that physiologically relevant concentrations of both Ca2+/Mg2+ interacted with the same locus. Moreover, the D67A/D73A N-terminal cTnC construct in which two coordinating residues within site II were removed was found to have significantly reduced affinity for both cations. In addition, 1 mM Mg2+ caused a 1.4-fold lower affinity for Ca2+. These experiments strongly suggest that cytosolic-free Mg2+ occupies a significant population of the available site II. Interaction of Mg2+ with site II of cTnC likely has important functional consequences for the heart both at baseline as well as in diseased states that decrease or increase the availability of Mg2+, such as secondary hyperparathyroidism or ischemia, respectively.

Nucleation plays a critical role in many natural and technological processes, and nucleation control requires detailed understanding of nucleation process at atomic level. In this study, we investigate the atomistic mechanism of heterogeneous nucleation in generic systems of liquid/substrate with positive lattice misfit (the solid has larger atomic spacing than the substrate) using molecular dynamics (MD) simulations. We found that heterogeneous nucleation process in such systems can be best described by a 3-layer nucleation mechanism: formation of the completely ordered first layer with an epitaxial relationship with the top surface of the substrate; formation of vacancies in the second layer to accommodate lattice misfit; and creation of a nearly perfect crystal plane of the solid in the third layer that demarcates the end of nucleation and the start of crystal growth. This 3-layer nucleation process creates a 2D nucleus (a plane of the solid phase), which contrasts with the hemisphere of the solid (a 3D nucleus) in the classical nucleation theory (CNT). It is expected that this 3-layer nucleation mechanism will provide new insight for nucleation control through effective manipulation of the liquid/substrate interface.

The pathogenesis of the SARS-CoV-2 virus initiates through recognition of the angiotensin-converting enzyme 2 (ACE2) receptor of the host cells by the receptor-binding domain (RBD) located at the spikes of the virus. Here, using molecular dynamics simulations, we have demonstrated the allosteric crosstalk within the RBD in the apo- and the ACE2 receptor-bound states, revealing the contribution of the dynamics-based correlated motions and the electrostatic energy perturbations to this crosstalk. While allostery, based on correlated motions, dominates inherent distal communication in the apo-RBD, the electrostatic energy perturbations determine favorable pairwise crosstalk within the RBD residues upon binding to ACE2. Interestingly, the allosteric path is composed of residues which are evolutionarily conserved within closely related coronaviruses, pointing toward the biological relevance of the communication and its potential as a target for drug development.

&beta, Glucosidases and &beta, mannosidases hydrolyze substrates that differ only in the epimer of the nonreducing terminal sugar moiety, but most such enzymes show a strong preference for one activity or the other. Rice Os3BGlu7 and Os7BGlu26 &beta, glycosidases show a less strong preference, but Os3BGlu7 and Os7BGlu26 prefer glucosides and mannosides, respectively. Previous studies of crystal structures with glucoimidazole (GIm) and mannoimidazole (MIm) complexes and metadynamic simulations suggested that Os7BGlu26 hydrolyzes mannosides via the B2,5 transition state (TS) conformation preferred for mannosides and glucosides via their preferred 4H3/4E TS conformation. However, MIm is weakly bound by both enzymes. In the present study, we found that MIm was not bound in the active site of crystallized Os3BGlu7, but GIm was tightly bound in the &minus, 1 subsite in a 4H3/4E conformation via hydrogen bonds with the surrounding residues. One-microsecond molecular dynamics simulations showed that GIm was stably bound in the Os3BGlu7 active site and the glycone-binding site with little distortion. In contrast, MIm initialized in the B2,5 conformation rapidly relaxed to a E3/4H3 conformation and moved out into a position in the entrance of the active site, where it bound more stably despite making fewer interactions. The lack of MIm binding in the glycone site in protein crystals and simulations implies that the energy required to distort MIm to the B2,5 conformation for optimal active site residue interactions is sufficient to offset the energy of those interactions in Os3BGlu7. This balance between distortion and binding energy may also provide a rationale for glucosidase versus mannosidase specificity in plant &beta, glycosidases.

There is a strong negative correlation between the polyglutamine (polyQ) domain length (Q-length) in the intrinsically disordered Huntingtin protein (Htt) exon-1 and the age of onset of Huntington’s disease (HD). PolyQ of Q-length longer than 40 has the propensity of forming very compact aggregate structures, leading to HD at full penetrance. Recent advances in nanobiotechnology has provided a new platform for the development of novel diagnosis and therapeutics. Here, we explore the possibility of utilizing 2D-nanomaterials to inhibit the formation of supercompact polyQ structures through the so-called “folding-upon-binding” where the protein structure is dependent on the binding substrate. Using molecular dynamics simulations, we characterize two polyQ peptides with Q-length of 22 (Q22, normal length) and 46 (Q46, typical length causing HD) binding to both graphene and molybdenum disulfide (MoS2) nanosheets, which have been applied as antibacterial or anticancer agents. Upon binding, Q22 unfolds and elongates on both grapheme and MoS2 surfaces, regardless of its initial conformation, with graphene showing slightly stronger effect. In contrast, initially collapsed Q46 remains mostly collapsed within our simulation time on both nanosheets even though they do provide some “stretching” to Q46 as well. Further analyses indicate that the hydrophobic nature of graphene/MoS2 promotes the stretching of polyQ on nanosheets. However, there is strong competition with the intra-polyQ interactions (mainly internal hydrogen bonds) leading to the disparate folding/binding behaviors of Q22 and Q46. Our results present distinct Q-length specific behavior of the polyQ domain upon binding to two types of 2D-nanomaterials which holds clinical relevance for Huntington's disease.

Protein design is able to create artificial proteins with advanced functions, and computer simulation plays a key role in guiding the rational design. In the absence of structural evidence for cytoglobin (Cgb) with an intramolecular disulfide bond, we recently designed a de novo disulfide bond in myoglobin (Mb) based on structural alignment (i.e., V21C/V66C Mb double mutant). To provide deep insight into the regulation role of the Cys21-Cys66 disulfide bond, we herein perform molecular dynamics (MD) simulation of the fluoride&ndash, protein complex by using a fluoride ion as a probe, which reveals detailed interactions of the fluoride ion in the heme distal pocket, involving both the distal His64 and water molecules. Moreover, we determined the kinetic parameters of fluoride binding to the double mutant. The results agree with the MD simulation and show that the formation of the Cys21-Cys66 disulfide bond facilitates both fluoride binding to and dissociating from the heme iron. Therefore, the combination of theoretical and experimental studies provides valuable information for understanding the structure and function of heme proteins, as regulated by a disulfide bond. This study is thus able to guide the rational design of artificial proteins with tunable functions in the future.

The inhibition of human angiotensin I converting enzyme (ACE) has been regarded as a promising approach for the treatment of hypertension. Despite research attempts over many years, our understanding the mechanisms of activation and inhibition of ACE is still far from complete. Here, we present results of all atom molecular dynamics simulations of ACE with and without ligands. Two types of inhibitors, competitive and mixed non-competitive, were used to model the ligand bound forms. In the absence of a ligand the simulation showed spontaneous large hinge-bending motions of multiple conversions between the closed and open states of ACE, while the ligand bound forms were stable in the closed state. Our simulation results imply that the equilibrium between pre-existing backbone conformations shifts in the presence of a ligand. The hinge-bending motion of ACE is considered as an essential to the enzyme function. A mechanistic model of activation and the inhibition may provide valuable information for novel inhibitors of ACE.

Recent studies have pointed the role of angiotensin-converting enzyme-II (ACE2) in mediating the entry of SARS-CoV-2 to the host cell by binding to the receptor-binding domain (RBD) of viral spike protein, and successive priming by cellular proteases initiates the infection. SARS-CoV replication rate and disease severity is controlled by the binding affinity of RBD with ACE2. To understand, how mutations in the conserved residues of RBD affect the molecular interaction with ACE2, we generated five alanine mutants i.e. Y449A, N487A, Y489A, N501A and Y505A in the receptor binding motif (RBM) of the ACE2-RBD SARS-CoV-2 complex (PDB: 6M0J). Computational site directed mutagenesis induced dynamics in wild-type and mutant complexes were extensively studied through all-atoms molecular dynamics (MD) simulations of 150 ns. In silico mutational analysis revealed loss of important intermolecular hydrogen bonds and other non-bonded contacts, critical for molecular recognition of SARS-CoV-2 RBD to ACE2, which is well supported by saturation mutagenesis study of binding interface residues. MD simulations results showed that RBM motif is flexible, where mutant residues are relatively more mobile than corresponding wild-type residues. Global motion analysis through principal component studies revealed that RBD exhibits protuberant in-ward motion towards the human ACE2 binding interface which may be crucial for molecular interaction. Conclusively, the present finding are in congruence with previous experimental reports and provides detailed information on the structural basis of receptor binding by human SARS-CoV-2, which will crucial for the development of novel inhibitors or drugs to combat against SARS-CoV-2. Communicated by Ramaswamy H. Sarma

We report data on the structural dynamics of the neuropeptide Y (NPY) G-protein-coupled receptor (GPCR) type 1 (Y1R), a typical representative of class A peptide ligand GPCRs, using a combination of solid-state NMR and molecular dynamics (MD) simulation. First, the equilibrium dynamics of Y1R were studied using 15N-NMR and quantitative determination of 1H-13C order parameters through the measurement of dipolar couplings in separated-local-field NMR experiments. Order parameters reporting the amplitudes of the molecular motions of the C-H bond vectors of Y1R in DMPC membranes are 0.57 for the C&alpha, sites and lower in the side chains (0.37 for the CH2 and 0.18 for the CH3 groups). Different NMR excitation schemes identify relatively rigid and also dynamic segments of the molecule. In monounsaturated membranes composed of longer lipid chains, Y1R is more rigid, attributed to a higher hydrophobic thickness of the lipid membrane. The presence of an antagonist or NPY has little influence on the amplitude of motions, whereas the addition of agonist and arrestin led to a pronounced rigidization. To investigate Y1R dynamics with site resolution, we conducted extensive all-atom MD simulations of the apo and antagonist-bound state. In each state, three replicas with a length of 20 &mu, s (with one exception, where the trajectory length was 10 &mu, s) were conducted. In these simulations, order parameters of each residue were determined and showed high values in the transmembrane helices, whereas the loops and termini exhibit much lower order. The extracellular helix segments undergo larger amplitude motions than their intracellular counterparts, whereas the opposite is observed for the loops, Helix 8, and termini. Only minor differences in order were observed between the apo and antagonist-bound state, whereas the time scale of the motions is shorter for the apo state. Although these relatively fast motions occurring with correlation times of ns up to a few &micro, s have no direct relevance for receptor activation, it is believed that they represent the prerequisite for larger conformational transitions in proteins.

We studied and compared in detail two published crystallographic structures of ADP/ATP carrier (AAC) protein in matrix open (m-state) and cytoplasmic open (c-state) state with and without cardiolipin (CDL) molecules embedded in the DOPC bilayer. We simulated and analyzed both states and obtained results show that both states are stable and water impermeable. In comparison with c-state, m-state is more occluded and has longer water impermeable area. Protein in both states, especially in the m- state, is more similar to crystallographic structure and less flexible in simulations when CDL molecules are present. This stabilization of m-state with CDL molecules sheds new light on AAC protein stability, representing a crucial factor for this structure which is significantly less stable compared with c- state.

Calcite and magnesite are important mineral constituents of the earth’s crust. In aqueous environments, these carbonates typically expose their most stable cleavage plane, the (10.4) surface. It is known that these surfaces interact with a large variety of organic molecules, which can result in surface restructuring. This process is decisive for the formation of biominerals. With the development of 3D atomic force microscopy (AFM) it is now possible to image solid–liquid interfaces with unprecedented molecular resolution. However, the majority of 3D AFM studies have been focused on the arrangement of water at carbonate surfaces. Here, we present an analysis of the assembly of ethanol – an organic molecule with a single hydroxy group – at the calcite and magnesite (10.4) surfaces by using high-resolution 3D AFM and molecular dynamics (MD) simulations. Within a single AFM data set we are able to resolve both the first laterally ordered solvation layer of ethanol on the calcite surface as well as the following solvation layers that show no lateral order. Our experimental results are in excellent agreement with MD simulations. The qualitative difference in the lateral order can be understood by the differing chemical environment: While the first layer adopts specific binding positions on the ionic carbonate surface, the second layer resides on top of the organic ethyl layer. A comparison of calcite and magnesite reveals a qualitatively similar ethanol arrangement on both carbonates, indicating the general nature of this finding.

Intrinsically Disordered Peptides and Proteins (IDPs) in solution can span a broad range of conformations that often are hard to characterize by both experimental and computational methods. However, obtaining a significant representation of the conformational space is important to understand mechanisms underlying protein functions such as partner recognition. In this work, we investigated the behavior of the Sic1 Kinase-Inhibitor Domain (KID) in solution by Molecular Dynamics (MD) simulations. Our results point out that application of common descriptors of molecular shape such as Solvent Accessible Surface (SAS) area can lead to misleading outcomes. Instead, more appropriate molecular descriptors can be used to define 3D structures. In particular, we exploited Weighted Holistic Invariant Molecular (WHIM) descriptors to get a coarse-grained but accurate definition of the variegated Sic1 KID conformational ensemble. We found that Sic1 is able to form a variable amount of folded structures even in absence of partners. Among them, there were some conformations very close to the structure that Sic1 is supposed to assume in the binding with its physiological complexes. Therefore, our results support the hypothesis that this protein relies on the conformational selection mechanism to recognize the correct molecular partners.

Here, we report on a computational comparison of the receptor-binding domains (RBDs) on the spike proteins of severe respiratory syndrome coronavirus‐2 (SARS-CoV-2) and SARS-CoV in free forms and as complexes with angiotensin-converting enzyme 2 (ACE2) as their receptor in humans. The impact of 42 mutations discovered so far on the structure and thermodynamics of SARS-CoV-2 RBD was also assessed. The binding affinity of SARS-CoV-2 RBD for ACE2 is higher than that of SARS-CoV RBD. The binding of COVA2-04 antibody to SARS-CoV-2 RBD is more energetically favorable than the binding of COVA2-39, but also less favorable than the formation of SARS-CoV-2 RBD-ACE2 complex. The net charge, the dipole moment and hydrophilicity of SARS-CoV-2 RBD are higher than those of SARS-CoV RBD, producing lower solvation and surface free energies and thus lower stability. The structure of SARS-CoV-2 RBD is also more flexible and more open, with a larger solvent-accessible surface area than that of SARS-CoV RBD. Single-point mutations have a dramatic effect on distribution of charges, most prominently at the site of substitution and its immediate vicinity. These charge alterations alter the free energy landscape, while X→F mutations exhibit a stabilizing effect on the RBD structure through π stacking. F456 and W436 emerge as two key residues governing the stability and affinity of the spike protein for its ACE2 receptor. These analyses of the structural differences and the impact of mutations on different viral strains and members of the coronavirus genera are an essential aid in the development of effective therapeutic strategies. Communicated by Ramaswamy H. Sarma, Graphical abstract

Tubitak Directorate of Research and Development Resources of Kadir Has University Histone deacetylase (HDAC) 10 has been implicated in the pathology of various cancers and neurodegenerative disorders, making the discovery of novel inhibitors of the isoform an important endeavor. However, the unavailability of crystallographic structure of human HDAC10 (hHDAC10) hinders structure-based drug design effort. Previously, we reported the homology modeled structure of human HDAC10 built using the crystallographic structure of Danio rerio (zebrafish) HDAC10 (zHDAC10) (Protein Data Bank (PDB) ID; 5TD7, released on 24 May 2017) as a template. Here, in continuation with our study, both hHDAC10 and zHDAC10, and their respective complexes with trichostatin A (TSA), quisinostat, and the native ligand (in 5TD7), 7-[(3-aminopropyl)amino]-1,1,1-trifluoroheptane-2,2-diol (PDB ID; FKS) were submitted to 100 ns-long unrestrained molecular dynamics (MD) simulations. Comparative analyses of the MD trajectories revealed that zHDAC10 and its complexes displayed higher stability than hHDAC10 and its corresponding complexes over time. Nonetheless, docking of active and inactive set molecules revealed that more reliable conformations of hHDAC10 could be obtained at an extended time period. This study may shed more light on the reliability of hHDAC10 modeled structure for use in selective inhibitor design. Communicated by Ramaswamy H. Sarma

The outbreak of a new coronavirus (SARS-CoV-2) was first identified in Wuhan, People's Republic of China, in 2019, which leads to a severe, life-threatening form of pneumonia (COVID-19). Research scientists all around the world have been trying to find small molecule drugs to treat COVID-19. In the present study, a conserved macrodomain, ADP Ribose phosphatase (ADRP), of a critical non-structural protein (Nsp3) in all coronaviruses was probed using large-scale Molecular Dynamics (MD) simulations to identify novel inhibitors. In our virtual screening workflow, the recently-solved X-ray complex structure, 6W6Y, in addition to a substrate-mimics was used to screen 17 million ZINC15 compounds using drug property filters and Glide docking scores. The top twenty output compounds each underwent 200 ns MD simulations (i.e. 20 × 200 ns) to validate their individual stability as potential inhibitors. Eight out of the twenty compounds showed stable binding modes in the MD simulations, as well as favorable drug properties from our predications. Therefore, our computational data suggest that the resulting top eight out of twenty compounds could potentially be novel inhibitors of ADRP for ADRP of SARS-CoV-2.

Graphical abstract Theoretical simulations reveal that phosphorene nanosheet has a considerable ability as a sensor for SARS-CoV-2 detection., Due to the dramatic increase in the number of patients with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), designing new selective and sensitive sensors for the detection of this virus is of importance. In this research, by employing full atomistic molecular dynamics (MD) simulations, the interactions of the receptor-binding domain (RBD) of the SARS-CoV-2 with phosphorene and graphene nanosheets were analyzed to investigate their sensing ability against this protein. Based on the obtained results, the RBD interactions with the surface of graphene and phosphorene nanosheets do not have important effects on the folding properties of the RBD but this protein has unique dynamical behavior against each nanostructure. In the presence of graphene and phosphorene, the RBD has lower stability because due to the strong interactions between RBD and these nanostructures. This protein spreads on the surface and has lower structural compaction, but in comparison with graphene, RBD shows greater stability on the surface of the phosphorene nanosheet. Moreover, RBD forms a more stable complex with phosphorene nanosheet in comparison with graphene due to greater electrostatic and van der Waals interactions. The calculated Gibbs binding energy for the RBD complexation process with phosphorene and graphene are −200.37 and −83.65 kcal mol−1, respectively confirming that phosphorene has higher affinity and sensitivity against this protein than graphene. Overall, the obtained results confirm that phosphorene can be a good candidate for designing new nanomaterials for selective detection of SARS-CoV-2.

The binding free-energy profile of methanol to Candida antarctica lipase B (CALB) was calculated at infinite dilution and at a finite methanol concentration of 6.1 M using umbrella sampling molecular dynamics simulations with the OPLS all-atom force field. An additional validation of the results was performed by employing alchemical double decoupling simulations. The binding free-energy profiles have been used in a related research article to validate free-energy profiles obtained from direct counting simulations with the aim to use the kinetic information encoded in the latter. The data provided in this work will be useful to study concentration effects on binding, to test alternative free energy methods or to use the proposed simulation protocol for related systems.