Comparative studies of structure and dynamics of caprine, leporine, ovine, and equine serum albumins.

Serum albumin (SA) is the most prevalent protein found in blood. Human albumin was used as an albumin substitute in hypoalbuminemia pets due to high sequence similarity. SAs from furry animals were also reported to be the major indoor allergens. Sensitizing to one of SAs coupled with high sequence identity can lead to cross-reactive antibodies in allergic individuals. Thus, understanding the structural and dynamic characters of SAs is crucial for not only albumin substitution but also allergen therapy. Herein, Molecular dynamics (MD) simulations were performed to elucidate the structural and dynamic dissimilarity and similarity of economic animals [equine (ESA), caprine (CASA), ovine (OSA), and leporine (LSA)] to albumins from human (HSA), bovine (BSA), porcine (PSA), and pets [cat (FSA) and dog (CSA)]. The aim is to explore the feasibility of HSA substitution and understand how albumins cause the cross-reactivity. Generally, all albumins studied here show the scissoring motion like other mammalian albumins. The uniqueness of each albumin is defined by different sequence identity of domain I. Also, the drug binding affinity of studied albumins differs from HSA, CSA, FSA, BSA, and PSA. Especially, LSA displays the most deviated behavior from the group. So, such albumin may not be suitable for albumin therapy for pets and humans. CASA, OSA, and ESA share similar characteristics, therefore it is possible to use them to monitor the osmotic pressure among their species, but the allergenic response must be seriously considered. An insight obtained here can be useful to develop albumin therapy and understand clinical allergy.Communicated by Ramaswamy H. Sarma.

Binding of Apo and Glycated Human Serum Albumins to an Albumin-Selective Aptamer-Bound Graphene Quantum Dot Complex.

Diabetes mellitus is a chronic metabolic disease involving continued elevated blood glucose levels. It is a leading cause of mortality and reduced life expectancy. Glycated human serum albumin (GHSA) has been reported to be a potential diabetes biomarker. A nanomaterial-based aptasensor is one of the effective techniques to detect GHSA. Graphene quantum dots (GQDs) have been widely used in aptasensors as an aptamer fluorescence quencher due to their high biocompatibility and sensitivity. GHSA-selective fluorescent aptamers are first quenched upon binding to GQDs. The presence of albumin targets results in the release of aptamers to albumin and consequently fluorescence recovery. To date, the molecular details on how GQDs interact with GHSA-selective aptamers and albumin remain limited, especially the interactions of an aptamer-bound GQD (GQDA) with an albumin. Thus, in this work, molecular dynamics simulations were used to reveal the binding mechanism of human serum albumin (HSA) and GHSA to GQDA. The results show the rapid and spontaneous assembly of albumin and GQDA. Multiple sites of albumins can accommodate both aptamers and GQDs. This suggests that the saturation of aptamers on GQDs is required for accurate albumin detection. Guanine and thymine are keys for albumin-aptamer clustering. GHSA gets denatured more than HSA. The presence of bound GQDA on GHSA widens the entrance of drug site I, resulting in the release of open-chain glucose. The insight obtained here will serve as a base for accurate GQD-based aptasensor design and development.

Exploring the structural and dynamic differences between human carnosinase I (CN1) and II (CN2).

Human carnosinases (CNs) are dimeric dipeptidases in the metallopeptidase M20 family. Two isoforms of carnosinases (Zn2+ -containing carnosinase 1 (CN1) found in serum and Mn2+ -carnosinase 2 (CN2) in tissue) were identified. Both CNs cleave histidine-containing (Xaa-His) dipeptides such as carnosine where CN2 was found to accept a broader spectrum of substrates. A loss of CN function, resulting in a high carnosine concentration, reduces risk for diabetes and neurological disorders. Although several studies on CN activities and its Michaelis complex were conducted, all shed the light on CN1 activity where the CN2 data is limited. Also, the molecular details on CN1 and CN2 similarity and dissimilarity in structure and function remain unclear. Thus, in this work, molecular dynamics (MD) simulations were employed to study structure and dynamics of human CN1 and CN2 in comparison. The results show that the different catalytic ability of both CNs is due to their pocket size and environment. CN2 can accept a wider range of substrate due to the wider mouth of a binding pocket. The L1 loop seems to play a role in gating activity. Comparing to CN2, CN1 provides more electronegative entrance, more wettability, and higher stability of catalytic metal ion-pair in the active site which allow more efficient water-mediated catalysis. The microscopic understanding obtained here can serve as a basis for CN inhibition strategies resulting in higher carnosine levels and consequently mitigating complications associated with diseases such as diabetes and neurological disorder.

What make malarial adenosine deaminase from PLASMODIUM VIVAX recognise adenosine and 5'-methylthioadenosine: simulation studies.

Malaria is a life-threatening disease in humans caused by Plasmodium parasites. Plasmodium vivax (P. vivax) is one of the prevalent species found worldwide. An increase in an anti-malarial drug resistance suggests the urgent need for new drugs. Zn2+-containing adenosine deaminase (ADA) is a promising drug target because the ADA inhibition is fatal to the parasite. Malarial ADA accepts both adenosine (ADN) and 5'-methylthioadenosine (MTA) as substrates. The understanding of the substrate binding becomes crucial for an anti-malarial drug development. In this work, ADA from P. vivax (pvADA) is of interest due to its prevalence worldwide. The binding of ADN and MTA are studied here using Molecular Dynamics (MD) simulations. Upon binding, the open and closed states of pvADA are captured. The displacement of α7, linking loops of ß3/α12, ß4/α13, ß5/α15, and α10/α11 is involved in the cavity closure and opening. Also, the inappropriate substrate orientation induces a failure in a complete cavity closure. Interactions with D46, D172, S280, D310, and D311 are important for ADN binding, whereas only hydrogen bonds with D172 and D311 are sufficient to anchor MTA inside the pocket. No Zn2+-coordinated histidine residues is acquired for substrate binding. D172 is found to play a role in ribose moiety recognition, while D311 is crucial for trapping the amine group of an adenine ring towards the Zn2+ site. Comparing between ADN and MTA, the additional interaction between D310 and an amine nitrogen on ADN supports a tighter fit that may facilitate the deamination.Communicated by Ramaswamy H. Sarma.

Exploring the binding modes of cordycepin to human adenosine deaminase 1 (ADA1) compared to adenosine and 2'-deoxyadenosine.

Cordycepin (3'-deoxyadenosine, abbreviated as COR) from Cordyceps shows a wide range of pharmacological activities, including antioxidant and anticancer effects, therefore representing a potential alternative medicine. However, COR has a short half-life in the human body, where it is metabolized by adenosine deaminase 1 (ADA1). ADA1 helps regulate adenosine levels by deaminating excess adenosine (ADE) and its derivatives, such as 2'-deoxyadenosine (DEO). Understanding binding mechanisms of ADA1 with COR in comparison with its other substrates will play a vital role in improving the bioactivity and lifetime of COR for commercial medicinal use. Recently, the first structure of human ADA1 in complex with DEO was solved. We therefore employed molecular dynamics (MD) simulations to predict structures and dynamics of ADA1 complexing with ADE, DEO, and COR in comparison to a ligand-free (LF) structure. Our data reveal that a large and highly water-exposed binding pocket of ADA1 is responsible for ligand translocation and reorientation. Two possible binding locations (site1 and site2) are identified. The binding affinities of the ligands are ADE > COR > DEO. Furthermore, the movements of two loop regions at the binding pocket entrance, residues 183-193 and 215-230, contribute to gating activity.

Dynamic and structural insights into tick serpin from Ixodes ricinus.

Ixodid ticks have a crucial impact on people and domestic animals worldwide. These parasites also pose a serious threat to livestock. To date, vaccination of hosts against ticks is a safer, more sustainable alternative to chemical control of ticks and the disease agents they transmit. Because of their roles in tick physiology, serpins (serine protease inhibitors) from tick saliva are among the candidates for anti-tick vaccines. Inhibitory serpins employ a suicide inhibition mechanism to inhibit proteases, where the serpin reactive centre loop (RCL) is cleaved, by the targeted protease, and then inserted into the main ß-sheet of the serpin. This causes a massive conformational change called the 'stressed to relaxed' (S→R) transition, leading to the breakdown of serpin into two regions (core domain and cleaved polypeptide). Recently, the first tick serpin crystal structure from Ixodes ricinus in R-state was reported. We thus employed molecular dynamics simulations to better understand serpin structure and dynamics in atomic detail. Overall, R-state serpin showed high rigidity, especially the core domain. The most flexible region is the terminal of the cleaved polypeptide, due to its high-water exposure, while the rest of the cleaved polypeptide is stably trapped behind the core domain. T363, D367 and N375 are found to play a vital role in protein-protein attachment. This finding can be used to explain the high stability of the R-state serpin at the atomic level and provides insight into this tick serpin which will be useful for rational anti-tick vaccine development. AbbreviationsMDMolecular DynamicsRCLReactive centre loopCommunicated by Ramaswamy H. Sarma.

Dynamics of human defensin 5 (HD5) self-assembly in solution: Molecular simulations/insights.

Human α -defensin 5 (HD5) is a 32-residue cysteine-rich host-defense peptide that exhibits broad-spectrum antimicrobial activity and plays an essential role in innate immunity in the human gut and other organ systems. Although its antimicrobial mechanism of action remains unclear, the high salt concentration seems to attenuate the antimicrobial function of HD5 via an unknown mechanism. In this work, we employ Molecular Dynamics (MD) simulations to analyse the oligomerization behaviour of HD5 when exposed to different salt concentration. We demonstrate that the presence of salt, such as sodium chloride (NaCl), promotes HD5 to form higher-order oligomers (up to heptamers) in our simulations. In addition, we also analyse the electrostatic interactions between the two Glu residues (E14 and E21) and their neighbouring residues. Our data confirm that the E14 residue is essential for the structural integrity, whereas the E21 residue contributes to the dimerization of HD5, suggesting that these Glu residues are important for the antimicrobial function of this peptide.

Probing the binding affinities of imipenem and ertapenem for outer membrane carboxylate channel D1 (OccD1) from P. aeruginosa: simulation studies.

Pseudomonas aeruginosa is an important nosocomial human pathogen. The major difficulty in the fight against this pathogen is the relative impermeability of its outer membrane (OM). Only specific substrates can penetrate through the OM of P. aeruginosa via substrate-specific porins, so this has become one of the most problematic drug-resistant pathogens. Carbapenems are the most effective drugs for treating P. aeruginosa infections. One such carbapenem that is applied in cases of P. aeruginosa infection is imipenem (IMI), which uses outer membrane carboxylate channel D1 (OccD1) as a point of entry into the pathogen. Unlike IMI, ertapenem (ERTA, another carbapenem) shows only weak activity towards P. aeruginosa, as it is blocked from penetrating through the OM. However, it is currently unclear as to why IMI is allowed to pass through the OM while ERTA is not. Therefore, we conducted molecular dynamics (MD) simulations to elucidate the behavior of these drugs inside OccD1 as compared to the ligand-free state. We discovered another possible binding site in the constriction region close to the side-pore opening. Both drugs employ the core lactam part to tether themselves to the binding site, whereas the tail governs the direction of permeation. L132 and F133 appear to be involved in interactions that are key to core attachment. At least four hydrogen bonds are required for drug binding. The direction of motion of L2 also plays a role: inward flipping traps IMI in the constriction area, while a shift of L2 towards the membrane brings ERTA into contact with more water, which prompts the expulsion of ERTA to the mouth of the channel protein. The opening of L2 seems to facilitate the rejection of ERTA.

Why do the outer membrane proteins OmpF from E. coli and OprP from P. aeruginosa prefer trimers? Simulation studies.

Porins are water-filled protein channels across the outer membrane of gram-negative bacteria. They facilitate the uptake of nutrients and essential ions. Solutes are filtered by a constriction loop L3 at the mid of a pore. Porins are heat-stable and resistant to toxic agents and detergents. Most porins are trimer, but no clear explanation why trimeric form is preferable. In this work, we thus studied effects of oligomerization on porin structure and function in microscopic detail. A well-studied OmpF (general porin from Escherichia coli) and well-characterised OprP (phosphate-specific pore from Pseudomonas aeruginosa) are used as samples from 2 types of porins found in gram-negative bacteria. MD simulations of trimeric and monomeric pores in pure water and 1M NaCl solution were performed. With a salt solution, the external electric field was applied to mimic a transmembrane potential. Expectedly, OprP is more stable than OmpF. Interestingly, being a monomer turns OmpF into an anion-selective pore. The dislocation of D113's side chain on L3 in OmpF causes the disruption of cation pathway resulting in the reduction of cation influx. In contrast, OprP's structure and function are less dependent on oligomeric states. Both monomeric and trimeric OprP can maintain their anion selectivity. Our findings suggest that trimerization is crucial for both structure and function of general porin OmpF, whereas being trimer in substrate-specific channel OprP supports a pore function.