Anisotropic remixing of a phase separated binary colloidal system with particles of different sizes in an external modulation.

We explore the phase behavior of a binary colloidal system under external spatially periodic modulation. We perform Monte Carlo simulations on a binary mixture of big and small repulsive Lennard-Jones particles with a diameter ratio of 2:1. We characterize structure by isotropic and anisotropic pair correlation functions, cluster size distribution, bond angle distribution, order parameter, and specific heat. We observe the demixing of the species in the absence of external modulation. However, the mixing of the species gets enhanced with increasing potential strength along with the alignment of the particles transverse to the modulation. The de-mixing order parameter shows discontinuity with increasing modulation strength, characterizing a first order phase transition. The peak in specific heat increases linearly with the size of the system. We also look into the dynamical behavior of the system via computing Mean Square Displacement (MSD) along both parallel and perpendicular directions to the modulation. We observe a decrease in the diffusion coefficient for both types of particles as we increase the strength of the modulation.

Fluctuation-Dominated Ligand Binding in Molten Globule Protein.

A molten globule (MG) state is an intermediate state of protein observed during the unfolding of the native structure. In MG states, milk protein α-lactalbumin (aLA) binds to oleic acid (OLA). This MG-aLA-OLA complex, popularly known as XAMLET, performs cytotoxic activities against cancer cell lines. However, the microscopic understanding of ligand recognition ability in the MG state of the protein has not yet been explored. Motivated by this, we explore the binding of bovine aLA with OLA using all-atom molecular dynamics (MD) simulations. We find the binding mode between MG-aLA and OLA using the conformational thermodynamics method. We also estimate the binding free energy using the umbrella sampling (US) method for both the MG state and the neutral state. We find that the binding free energy obtained from US is comparable with earlier experimental results. We characterize the dihedral fluctuations as the ligand is liberated from the active site of the protein using steered MD. The low energy fluctuations occur near the ligand binding site, which eventually transfer toward the Ca2+-binding site as the ligand is taken away from the protein.

Surface specific adsorption of glucose to ZnO.

ZnO is bio-safe and hence, may be a potential candidate for direct use as a glucose sensor. This requires an understanding of the interaction of glucose with four common surfaces, (101̄0), (112̄0), (0001) and (0001̄) of ZnO. We carry out molecular dynamics (MD) simulations enhanced by umbrella sampling of a glucose molecule in a solvent over a hydrated ZnO slab. The slab is obtained by quantum mechanical optimization. We observe that hydration layers formed above the surfaces affect the approach of glucose to the surfaces. Potential of mean force (PMF) calculations show that the (101̄0) surface shows the strongest adsorption of adsorption free energy -6.81 kJ mol-1 towards glucose. Thus, we offer a theoretical understanding on the interactions at the nano-bio junction of glucose and ZnO surfaces. Our study suggests that the (101̄0) surface may be used to fabricate a direct glucose sensor.

Effect on the conformations of the spike protein of SARS-CoV-2 due to mutation.

The spike protein of SARS-CoV-2 mediates receptor binding and cell entry and is the key immunogenic target for virus neutralization and the present attention of many vaccine layouts. It exhibits significant conformational flexibility. We study the structural fluctuations of spike protein among the most common mutations that appeared in the variant of concerns (VOC). We report the thermodynamics of conformational changes in mutant spike protein with respect to the wild-type from the distributions of the dihedral angles obtained from the equilibrium configurations generated via all-atom molecular dynamics simulations. We find that the mutation causes the increase in distance between the N-terminal domain and receptor binding domain, leading to an obtuse angle cosine θ distribution in the trimeric structure in spike protein. Thus, an increase in open state is conferred to the more infectious variants of SARS-CoV-2. The thermodynamically destabilized and disordered residues of receptor binding motif among the mutant variants of spike protein are proposed to serve as better binding sites for the host factor. We identify a short stretch of region connecting the N-terminal domain and receptor binding domain forming a linker loop where many residues undergo stabilization in the open state compared to the closed one.

Effect of biphosphate salt on dipalmitoylphosphatidylcholine bilayer deformation by Tat polypeptide.

Translocation of positively charged cell penetrating peptides (CPP) through cell membrane is important in drug delivery. Here we report all-atom molecular dynamics simulations to investigate how a biphosphate salt in a solvent affects the interaction of a CPP, HIV-1 Tat peptide with model dipalmitoylphosphatidylcholine (DPPC) lipid bilayer. Tat peptide has a large number of basic arginines and a couple of polar glutamines. We observe that in absence of salt, the basic residues of the polypeptide get localized in the vicinity of the membrane without altering the bilayer properties much; polypeptide induce local thinning of the bilayer membrane at the area of localization. In presence of biphosphate salt, the basic residues, dressed by the biphosphate ions, are repelled by the phosphate head groups of the lipid molecules. However, polar glutamine prefers to stay in the vicinity of the bilayer. This leads to larger local bilayer thickness at the contact point by the polar residue and non-uniform bilayer thickness profile. The thickness deformation of bilayer structure disappears upon mutating the polar residue, suggesting importance of the polar residue in bilayer deformation. Our studies point to control bilayer deformation by appropriate peptide sequence and solvent conditions.

Effective single component description of steady state structures of passive particles in an active bath.

We model a binary mixture of passive and active Brownian particles in two dimensions using the effective interaction between passive particles in the active bath. The activity of active particles and the size ratio of two types of particles are the two control parameters in the system. The effective interaction is calculated from the average force on two particles generated by the active particles. The effective interaction can be attractive or repulsive, depending on the system parameters. The passive particles form four distinct structural orders for different system parameters, viz., homogeneous structures, disordered cluster, ordered cluster, and crystalline structure. The change in structure is dictated by the change in nature of the effective interaction. We further confirm the four structures using a full microscopic simulation of active and passive mixture. Our study is useful to understand the different collective behavior in non-equilibrium systems.

Microscopic insight to specificity of metal ion cofactor in DNA cleavage by restriction endonuclease EcoRV.

Restriction endonucleases protect bacterial cells against bacteriophage infection by cleaving the incoming foreign DNA into fragments. In presence of Mg2+ ions, EcoRV is able to cleave the DNA but not in presence of Ca2+ , although the protein binds to DNA in presence of both metal ions. We make an attempt to understand this difference using conformational thermodynamics. We calculate the changes in conformational free energy and entropy of conformational degrees of freedom, like DNA base pair steps and dihedral angles of protein residues in Mg2+ (A)-EcoRV-DNA complex compared to Ca2+ (S)-EcoRV-DNA complex using all-atom molecular dynamics (MD) trajectories of the complexes. We find that despite conformational stability and order in both complexes, the individual degrees of freedom behave differently in the presence of two different metal ions. The base pairs in cleavage region are highly disordered in Ca2+ (S)-EcoRV-DNA compared to Mg2+ (A)-EcoRV-DNA. One of the acidic residues ASP90, coordinating to the metal ion in the vicinity of the cleavage site, is conformationally destabilized and disordered, while basic residue LYS92 gets conformational stability and order in Ca2+ (S) bound complex than in Mg2+ (A) bound complex. The enhanced fluctuations hinder placement of the metal ion in the vicinity of the scissile phosphate of DNA. Similar loss of conformational stability and order in the cleavage region is observed by the replacement of the metal ion. Considering the placement of the metal ion near scissile phosphate as requirement for cleavage action, our results suggest that the changes in conformational stability and order of the base pair steps and the protein residues lead to cofactor sensitivity of the enzyme. Our method based on fluctuations of microscopic conformational variables can be applied to understand enzyme activities in other protein-DNA systems.

Self-assembly in amphiphilic macromolecules with solvent exposed hydrophobic moieties.

Self-assembly by amphiphilic molecules with solvent exposed hydrophobic groups are relevant in biomolecular systems as well as in technological applications. Here we study such self-assembly in these systems using a model system of spherical particles having charge at core but solvent repelling surface, using Monte-Carlo simulations and mean field treatment. We find that solvophobicity mediated attraction leads aggregation, while electrostatic repulsions control stability of finite clusters. The aggregation threshold relates the parameters of two interactions through an algebraic dependence. The study also qualitatively explains experimental observations on aggregation of misfolded proteins and can be useful guide to tune stability of nm sized self-assembly in systems with exposed hydrophobic groups.

Molecular dynamics simulations on interaction between bacterial proteins: Implication on pathogenic activities.

We perform molecular dynamics simulation studies on interaction between bacterial proteins: an outer-membrane protein STY3179 and a yfdX protein STY3178 of Salmonella Typhi. STY3179 has been found to be involved in bacterial adhesion and invasion. STY3178 is recently biophysically characterized. It is a soluble protein having antibiotic binding and chaperon activity capabilities. These two proteins co-occur and are from neighboring gene in Salmonella Typhi-occurrence of homologs of both STY3178 and STY3179 are identified in many Gram-negative bacteria. We show using homology modeling, docking followed by molecular dynamics simulation that they can form a stable complex. STY3178 belongs to aqueous phase, while the beta barrel portion of STY3179 remains buried in DPPC bilayer with extra-cellular loops exposed to water. To understand the molecular basis of interaction between STY3178 and STY3179, we compute the conformational thermodynamics which indicate that these two proteins interact through polar and acidic residues belonging to their interfacial region. Conformational thermodynamics results further reveal instability of certain residues in extra-cellular loops of STY3179 upon complexation with STY3178 which is an indication for binding with host cell protein laminin.

Quantum chemical studies on anion specificity of CαNN motif in functional proteins.

Anion binding CαNN motif is found in functionally important regions of protein structures. This motif based only on backbone atoms from three adjacent residues, recognizes free sulphate or phosphate ion as well as phosphate groups in nucleotides and in a variety of cofactors. The mode of anion recognition and microscopic picture of binding interaction remains unclear. Here we perform self-consistent quantum chemical calculations considering sulphate and phosphate bound CαNN motif fragments from crystal structures of functional proteins in order to figure out microscopic basis of anion recognition. Our calculations indicate that stability and preference of the anion in the motif depends on the sequence of the motif. The stabilization energy is larger in case of polar residue containing motif fragment. Nitrogen atom of the polar residue of motif mainly participates in the coordination at the lowest energy levels. Anion replacement decreases stabilization energy along with coordination between motif atoms and oxygen atoms of anion shifted to higher energies, suggesting preference of the motif residues to specific anion. Our analysis may be helpful to understand microscopic basis of interaction between proteins and ionic species.

Anion induced conformational preference of Cα NN motif residues in functional proteins.

Among different ligand binding motifs, anion binding Cα NN motif consisting of peptide backbone atoms of three consecutive residues are observed to be important for recognition of free anions, like sulphate or biphosphate and participate in different key functions. Here we study the interaction of sulphate and biphosphate with Cα NN motif present in different proteins. Instead of total protein, a peptide fragment has been studied keeping Cα NN motif flanked in between other residues. We use classical force field based molecular dynamics simulations to understand the stability of this motif. Our data indicate fluctuations in conformational preferences of the motif residues in absence of the anion. The anion gives stability to one of these conformations. However, the anion induced conformational preferences are highly sequence dependent and specific to the type of anion. In particular, the polar residues are more favourable compared to the other residues for recognising the anion.

Hybrid simulation approach incorporating microscopic interaction along with rigid body degrees of freedom for stacking between base pairs.

Stacking interaction between the aromatic heterocyclic bases plays an important role in the double helical structures of nucleic acids. Considering the base as rigid body, there are total of 18 degrees of freedom of a dinucleotide step. Some of these parameters show sequence preferences, indicating that the detailed atomic interactions are important in the stacking. Large variants of non-canonical base pairs have been seen in the crystallographic structures of RNA. However, their stacking preferences are not thoroughly deciphered yet from experimental results. The current theoretical approaches use either the rigid body degrees of freedom where the atomic information are lost or computationally expensive all atom simulations. We have used a hybrid simulation approach incorporating Monte-Carlo Metropolis sampling in the hyperspace of 18 stacking parameters where the interaction energies using AMBER-parm99bsc0 and CHARMM-36 force-fields were calculated from atomic positions. We have also performed stacking energy calculations for structures from Monte-Carlo ensemble by Dispersion corrected density functional theory. The available experimental data with Watson-Crick base pairs are compared to establish the validity of the method. Stacking interaction involving A:U and G:C base pairs with non-canonical G:U base pairs also were calculated and showed that these structures were also sequence dependent. This approach could be useful to generate multiscale modeling of nucleic acids in terms of coarse-grained parameters where the atomic interactions are preserved. This method would also be useful to predict structure and dynamics of different base pair steps containing non Watson-Crick base pairs, as found often in the non-coding RNA structures. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 212-226, 2016.

Role of indirect readout mechanism in TATA box binding protein-DNA interaction.

Gene expression generally initiates from recognition of TATA-box binding protein (TBP) to the minor groove of DNA of TATA box sequence where the DNA structure is significantly different from B-DNA. We have carried out molecular dynamics simulation studies of TBP-DNA system to understand how the DNA structure alters for efficient binding. We observed rigid nature of the protein while the DNA of TATA box sequence has an inherent flexibility in terms of bending and minor groove widening. The bending analysis of the free DNA and the TBP bound DNA systems indicate presence of some similar structures. Principal coordinate ordination analysis also indicates some structural features of the protein bound and free DNA are similar. Thus we suggest that the DNA of TATA box sequence regularly oscillates between several alternate structures and the one suitable for TBP binding is induced further by the protein for proper complex formation.

In-silicon studies on hydration in EcoRI-cognate DNA complex.

Restriction endonucleases (REs) cleave DNA at specific site in presence of Mg2+ ion. Experiments further emphasize the role of hydration in metal ion specificity and sequence specificity of DNA cleavage. However, the relation between hydration and specificity has not been understood till date. This leads us to study via all-atom molecular dynamics (MD) simulations how the hydration around the scissile phosphate group changes in presence of Mg2+ and Ca2+ and depend on the DNA sequence. We observe the least number of hydrogen bonds around the scissile phosphate group in presence of Mg2+ ion. We further find that the hydrogen bonds decrease at the scissile phosphate on mutating one base pair in the cleavage region of the DNA in Mg2+ loaded EcoRI-DNA complex. We also perform steered MD simulations and observe that the rate of decrease of fraction of hydrogen bonds is slower in the mutated complex than the unmutated complex.

Immune escape facilitation by mutations of epitope residues in RdRp of SARS-CoV-2.

Mutations drive viral evolution and genome variability that causes viruses to escape host immunity and to develop drug resistance. SARS-CoV-2 has considerably higher mutation rate. SARS-CoV-2 possesses a RNA dependent RNA polymerase (RdRp) which helps to replicate its genome. The mutation P323L in RdRp is associated with the loss of a particular epitope (321-327) from this protein. We consider the effects of mutations in some of the epitope region including the naturally occurring mutation P323L on the structure of the epitope and their interface with paratope using all-atom molecular dynamics (MD) simulation studies. We observe that the mutations cause conformational changes in the epitope region by opening up the region associated with increase in the radius of gyration and intramolecular hydrogen bonds, making the region less accessible. Moreover, we study the conformational stability of the epitope region and epitope:paratope interface under the mutation from the fluctuations in the dihedral angles. We observe that the mutation renders the epitope and the epitope:paratope interface unstable compared to the corresponding wild type ones. Thus, the mutations may help in escaping antibody mediated immunity of the hostCommunicated by Ramaswamy H. Sarma.

Conformational stability and order of Hoogsteen base pair induced by protein binding.

Several experimental studies have shown that Hoogsteen (HG) base pair (bp) stabilizes in the presence of proteins. The molecular mechanism underlying this stabilization is not well known. This leads us to examine the stability of the HG bp in duplex DNA using all-atom molecular dynamics simulation in both the absence and presence of proteins. We use conformational thermodynamics to investigate the stability of a HG bp in duplex DNA at the molecular level. We compute the changes in the conformational free energy and entropy of DNA when DNA adopts a HG bp in its bp sequence rather than a Watson-Crick (WC) bp in both naked DNA and protein-bound DNA complex. We observe that the presence of proteins stabilizes and organizes the HG bp and the entire DNA duplex. Sugar-phosphate, sugar-base, and sugar-pucker torsion angles play key roles in stabilizing and ordering the HG bp in the protein-bound DNA complex.

Microscopic model on indoor propagation of respiratory droplets.

Indoor propagation of airborne diseases is yet poorly understood. Here, we theoretically study a microscopic model based on the motions of virus particles in a respiratory microdroplet, responsible for airborne transmission of diseases, to understand their indoor propagation. The virus particles are driven by a driving force that mimics force due to gushing of air by devices like indoor air conditioning along with the gravity. A viral particle within the droplet experiences viscous drag due to the droplet medium, force due to interfacial tension at the droplet boundary, the thermal forces and mutual interaction forces with the other viral particles. We use Brownian Dynamics (BD) simulations and scaling arguments to study the motion of the droplet, given by that of the center of mass of the viral assembly. The BD simulations show that in presence of the gravity force alone, the time the droplet takes to reach the ground level, defined by the gravitational potential energy being zero, from a vertical height H,tf∼γ-0.1 dependence, where γ is the interfacial tension. In presence of the driving force of magnitude F0 and duration τ0, the horizontal propagation length, Ymax from the source increase linearly with τ0, where the slope is steeper for larger F0. Our scaling analysis explains qualitatively well the simulation observations and show long-distance transmission of airborne respiratory droplets in the indoor conditions due to F0 ∼ nano-dyne.

Unraveling DPP4 Receptor Interactions with SARS-CoV-2 Variants and MERS-CoV: Insights into Pulmonary Disorders via Immunoinformatics and Molecular Dynamics.

Human coronaviruses like MERS CoV are known to utilize dipeptidyl peptidase 4 (DPP4), apart from angiotensin-converting enzyme 2(ACE2) as a potential co-receptor for viral cell entry. DPP4, the ubiquitous membrane-bound aminopeptidase, is closely associated with elevation of disease severity in comorbidities. In SARS-CoV-2, there is inadequate evidence for combination of spike protein variants with DPP4, and underlying adversity in COVID-19. To elucidate this mechanistic basis, we have investigated interaction of spike protein variants with DPP4 through molecular docking and simulation studies. The possible binding interactions between the receptor binding domain (RBD) of different spike variants of SARS-CoV-2 and DPP4 have been compared with interactions observed in the experimentally determined structure of the complex of MERS-CoV with DPP4. Comparative binding affinity confers that Delta-CoV-2: DPP4 shows close proximity with MERS-CoV:DPP4, as depicted from accessible surface area, radius of gyration and number of hydrogen bonding in the interface. Mutations in the delta variant, L452R and T478K directly participate in DPP4 interaction, enhancing DPP4 binding. E484K in alpha and gamma variants of spike protein is also found to interact with DPP4. Hence, DPP4 interaction with spike protein becomes more suitable due to mutation, especially due to L452R, T478K and E484K. Furthermore, perturbation in the nearby residues Y495, Q474 and Y489 is evident due to L452R, T478K and E484K, respectively. Virulent strains of spike protein are more susceptible to DPP4 interaction and are prone to be victimized in patients due to comorbidities. Our results will aid the rational optimization of DPP4 as a potential therapeutic target to manage COVID-19 disease severity.

A Mechanoelastic Glimpse on Hyaluronan-Coated Extracellular Vesicles.

Cancer cells secrete extracellular vesicles (EVs) covered with a carbohydrate polymer, hyaluronan (HA), linked to tumor malignancy. Herein, we have unravelled the contour lengths of HA on a single cancer cell-derived EV surface using single-molecule force spectroscopy (SMFS), which divulges the presence of low molecular weight HA (LMW-HA < 200 kDa). We also discovered that these LMW-HA-EVs are significantly more elastic than the normal cell-derived EVs. This intrinsic elasticity of cancer EVs could be directly allied to the LMW-HA abundance and associated labile water network on EV surface as revealed by correlative SMFS, hydration dynamics with fluorescence spectroscopy, and molecular dynamics simulations. This method emerges as a molecular biosensor of the cancer microenvironment.

2D-Molybdenum Disulfide-Derived Ion Source for Mass Spectrometry.

Generation of current or potential at nanostructures using appropriate stimuli is one of the futuristic methods of energy generation. We developed an ambient soft ionization method for mass spectrometry using 2D-MoS2, termed streaming ionization, which eliminates the use of traditional energy sources needed for ion formation. The ionic dissociation-induced electrokinetic effect at the liquid-solid interface is the reason for energy generation. We report the highest figure of merit of current generation of 1.3 A/m2 by flowing protic solvents at 22 µL/min over a 1 × 1 mm2 surface coated with 2D-MoS2, which is adequate to produce continuous ionization of an array of analytes, making mass spectrometry possible. Weakly bound ion clusters and uric acid in urine have been detected. Further, the methodology was used as a self-energized breath alcohol sensor capable of detecting 3% alcohol in the breath.