Mutations in exons 3 and 7 resulting in truncated expression of human ATP6V1B1 gene showing structural variations contributing to poor substrate binding-causative reason for distal renal tubular acidosis with sensorineural deafness

Distal renal tubular acidosis (dRTA) is an autosomal recessive syndrome results defect in either proximal tubule bicarbonate reabsorption or in distal tubule H⁺ secretion and is characterized by severe hyperchloraemic metabolic acidosis in childhood. dRTA is associated with functional variations in the ATP6V1B1 gene encoding β1 subunit of H⁺-ATPase, key membrane transporters for net acid excretion of α-intercalated cells of medullary collecting ducts. In the present study, a 13-year-old male patient suffering with nephropathy and sensorineural deafness was reported in the Department of Nephrology. We predicted improper functioning of ATP6V1B1 gene could be the reason for diseased condition. Therefore, exons 3, 4, and 7 contributing active site of ATP6V1B1 gene was amplified and sequenced (Accession numbers: KF571726, KM222653). The obtained sequences were BLAST searched against the wild type ATP6V1B1 gene which showed novel mutations c.307 A > G, c.308 C > A, c.310 C > G, c.704 T > C, c.705 G > T, c.709 A > G, c.710 A > G, c.714 G > A, c.716 C > A, c.717delC, c.722 C > G, c.728insG, c.741insT, c.753G > C. These mutations resulted in the expression of truncated protein terminating at Lys 209. The mutated ATP6V1B1structure superimposed with wild type showed extensive variations with RMSD 1.336 Å and could not bind to substrate ADP leading to non-functional ATPase. These results conclusively explain these mutations in ATP6V1B1 gene resulted in structural changes causing accumulation of H⁺ ions contributing to dRTA with sensorineural deafness.     

Solution structure and dynamics of human metallothionein-3 (MT-3)

Alzheimer's disease is characterized by progressive loss of neurons accompanied by the formation of intraneural neurofibrillary tangles and extracellular amyloid plaques. Human neuronal growth inhibitory factor, classified as metallothionein-3 (MT-3), was found to be related to the neurotrophic activity promoting cortical neuron survival and dendrite outgrowth in the cell culture studies. We have determined the solution structure of the alpha-domain of human MT-3 (residues 32-68) by multinuclear and multidimensional NMR spectroscopy in combination with the molecular dynamic simulated annealing approach. The human MT-3 shows two metal-thiolate clusters, one in the N-terminus (beta-domain) and one in the C-terminus (alpha-domain). The overall fold of the alpha-domain is similar to that of mouse MT-3. However, human MT-3 has a longer loop in the acidic hexapeptide insertion than that of mouse MT-3. Surprisingly, the backbone dynamics of the protein revealed that the beta-domain exhibits similar internal motion to the alpha-domain, although the N-terminal residues are more flexible. Our results may provide useful information for understanding the structure-function relationship of human MT-3.     

Comparative analysis of inner cavities and ligand migration in non-symbiotic AHb1 and AHb2

This study reports a comparative analysis of the topological properties of inner cavities and the intrinsic dynamics of non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana. The two proteins belong to the 3/3 globin fold and have a sequence identity of about 60%. However, it is widely assumed that they have distinct physiological roles. In order to investigate the structure–function relationships in these proteins, we have examined the bis-histidyl and ligand-bound hexacoordinated states by atomistic simulations using in silico structural models. The results allow us to identify two main pathways to the distal cavity in the bis-histidyl hexacoordinated proteins. Nevertheless, a larger accessibility to small gaseous molecules is found in AHb2. This effect can be attributed to three factors: the mutation Leu35(AHb1) → Phe32(AHb2), the enhanced flexibility of helix B, and the more favorable energetic profile for ligand migration to the distal cavity. The net effect of these factors would be to facilitate the access of ligands, thus compensating the preference for the fully hexacoordination of AHb2, in contrast to the equilibrium between hexa- and pentacoordinated species in AHb1. On the other hand, binding of the exogenous ligand introduces distinct structural changes in the two proteins. A well-defined tunnel is formed in AHb1, which might be relevant to accomplish the proposed NO detoxification reaction. In contrast, no similar tunnel is found in AHb2, which can be ascribed to the reduced flexibility of helix E imposed by the larger number of salt bridges compared to AHb1. This feature would thus support the storage and transport functions proposed for AHb2. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

Dimer Formation of a Stabilized Gβ1 Variant: A Structural and Energetic Analysis

In previous work, a strongly stabilized variant of the β1 domain of streptococcal protein G (Gβ1) was obtained by an in vitro selection method. This variant, termed Gβ1-M2, contains the four substitutions E15V, T16L, T18I, and N37L. Here we elucidated the molecular basis of the observed strong stabilizations. The contributions of these four residues were analyzed individually and in various combinations, additional selections with focused Gβ1 gene libraries were performed, and the crystal structure of Gβ1-M2 was determined. All single substitutions (E15V, T16L, T18I, and N37L) stabilize wild-type Gβ1 by contributions of between 1.6 and 6.0 kJ mol^− 1 (at 70 °C). Hydrophobic residues at positions 16 and 37 provide the major contribution to stabilization by enlarging the hydrophobic core of Gβ1. They also increase the tendency to form dimers, as shown by dependence on the concentration of apparent molecular mass in analytical ultracentrifugation, by concentration-dependent stability, and by a strongly increased van't Hoff enthalpy of unfolding. The 0.88-Å crystal structure of Gβ1-M2 and NMR measurements in solution provide the explanation for the observed dimer formation. It involves a head-to-head arrangement of two Gβ1-M2 molecules via six intermolecular hydrogen bonds between the two β strands 2 and 2′ and an adjacent self-complementary hydrophobic surface area, which is created by the T16L and N37L substitutions and a large 120° rotation of the Tyr33 side chain. This removal of hydrophilic groups and the malleability of the created hydrophobic surface provide the basis for the dimer formation of stabilized Gβ1 variants.     

Simulating structural and thermodynamic properties of carcinogen-damaged DNA

A pair of stereoisomeric covalent adducts to guanine in double-stranded DNA, derived from the reaction of mutagenic and tumorigenic metabolites of benzo[a]pyrene, have been well characterized structurally and thermodynamically. Both high-resolution NMR solution structures and an array of thermodynamic data are available for these 10S (+)- and 10R (-)-trans-anti -[BP]-N(2)-dG adducts in double-stranded deoxyoligonucleotides. The availability of experimentally well-characterized duplexes containing these two stereoisomeric guanine adducts provides an opportunity for evaluating the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method for computing thermodynamic properties from molecular dynamics ensembles. We have carried out 3-ns molecular dynamics simulations, using NMR solution structures as the starting models for the 10S (+)- and 10R (-)-trans-anti-dG adducts in a DNA duplex 11-mer using AMBER 6.0. We employed the MM-PBSA method to compute the free energies, enthalpies, and entropies of the two adducts. Our complete thermodynamic analysis agrees quite well with the full experimental thermodynamic characterization of these adducts, showing essentially equal stabilities of the two adducts. We also calculated the nuclear Overhauser effect (NOE) distances from the molecular dynamics trajectories, and compared them against the experimental NMR-derived NOE distances. Our results showed that the simulated structures are in good agreement with the NMR experimental NOE data. Furthermore, the molecular dynamics simulations provided new structural and biological insights. Specifically, the puzzling observation that the BP aromatic ring system in the 10S (+)-trans-anti-dG adduct is more exposed to the aqueous solvent than the 10R (-)-trans-anti-dG adduct, is rationalized in terms of the adduct structures. The structural and thermodynamic features of these stereoisomeric adducts are also discussed in relation to their reported low susceptibilities to nucleotide excision repair.     

Identification of a collapsed intermediate with non-native long-range interactions on the folding pathway of a pair of Fyn SH3 domain mutants by NMR relaxation dispersion spectroscopy

Recent 15N and 13C spin-relaxation dispersion studies of fast-folding mutants of the Fyn SH3 domain have established that folding proceeds through a low-populated on-pathway intermediate (I) where the central beta-sheet is at least partially formed, but without interactions between the NH2- and COOH-terminal beta-strands that exist in the folded state (F). Initial studies focused on mutants where Gly48 is replaced; in an effort to establish whether this intermediate is a general feature of Fyn SH3 folding a series of 15N relaxation experiments monitoring the folding of Fyn SH3 mutants N53P/V55L and A39V/N53P/V55L are reported here. For these mutants as well, folding proceeds through an on-pathway intermediate with similar features to those observed for G48M and G48V Fyn SH3 domains. However, the 15N chemical shifts extracted for the intermediate indicate pronounced non-native contacts between the NH2 and COOH-terminal regions not observed previously. The kinetic parameters extracted for the folding of A39V/N53P/V55L Fyn SH3 from the three-state folding model F<-->I<-->U are in good agreement with folding and unfolding rates extrapolated to zero denaturant obtained from stopped-flow experiments analyzed in terms of a simplified two-state folding reaction. The folding of the triple mutant was studied over a wide range of temperatures, establishing that there is no difference in heat capacities between F and I states. This confirms a compact folding intermediate structure, which is supported by the 15N chemical shifts of the I state extracted from the dispersion data. The temperature-dependent relaxation data simplifies data analysis because at low temperatures (< 25 degrees C) the unfolded state (U) is negligibly populated relative to I and F. A comparison between parameters extracted at low temperatures where the F<-->I exchange model is appropriate with those from the more complex, three-state model at higher temperatures has been used to validate the protocol for analysis of three-site exchange relaxation data.     

Solution structure of DNA/RNA hybrid duplex with C8-propynyl 2′-deoxyadenosine modifications: Implication of RNase H and DNA/RNA duplex interaction

Solution structures of DNA/RNA hybrid duplexes, d(GCGCA*AA*ACGCG): r(cgcguuuugcg)d(C) (designated PP57), containing two C8-propynyl 2′-deoxyadenosines (A*) and unmodified hybrid (designated U4A4) are solved. The C8-propynyl groups on 2′-deoxyadenosine perturb the local structure of the hybrid duplex, but overall the structure is similar to that of canonical DNA/RNA hybrid duplex except that Hoogsteen hydrogen bondings between A* and U result in lower thermal stability. RNase H is known to cleave RNA only in DNA/RNA hybrid duplexes. Minor groove widths of hybrid duplexes, sugar puckerings of DNA are reported to be responsible for RNase H mediated cleavage, but structural requirements for RNase H mediated cleavage still remain elusive. Despite the presence of bulky propynyl groups of PP57 in the minor groove and greater flexibility, the PP57 is an RNase H substrate. To provide an insight on the interactions between RNase H and substrates we have modeled Bacillus halodurans RNase H-PP57 complex, our NMR structure and modeling study suggest that the residue Gly(15) and Asn(16) of the loop residues between first β sheet and second β sheet of RNase HI of Escherichia coli might participate in substrate binding.     

Exploring the antimalarial potential of the methoxy-thiazinoquinone scaffold: Identification of a new lead candidate

A small library of antiplasmodial methoxy-thiazinoquinones, rationally designed on the model of the previously identified hit 1, has been prepared by a simple and inexpensive procedure. The synthetic derivatives have been subjected to in vitro pharmacological screening, including antiplasmodial and toxicity assays. These studies afforded a new lead candidate, compound 9, endowed with higher antiplasmodial potency compared to 1, a good selectivity index when tested against a panel of mammalian cells, no toxicity against RBCs, a synergistic antiplasmodial action in combination with dihydroartemisinin, and a promising inhibitory activity on stage V gametocyte growth. Computational studies provided useful insights into the structural requirements needed for the antiplasmodial activity of thiazinoquinone compounds and on their putative mechanism of action.     

Water transport in AQP0 aquaporin: molecular dynamics studies

A double lipid bilayer structure containing opposing tetramers of AQP0 aquaporin, in contact through extracellular face loop regions, was recently modeled using an intermediate-resolution map obtained by electron crystallographic methods. The pores of these water channels were found to be critically narrow in three regions and subsequently interpreted to be those of a closed state of the channel. The subsequent determination of a high-resolution AQP0 tetramer structure by X-ray crystallographic methods yielded a pore model featuring two of the three constrictions as noted in the EM work and water molecules within the channel pore. The extracellular-side constriction region of this AQP0 structure was significantly larger than that of the EM-based model and similar to that of the highly water permeable AQP1. The X-ray-based study of AQP0 however could not ascertain if the water molecules found in the pore were the result of water entering from one or both ends of the channel, nor whether water could freely pass through all constriction points. Additionally, this X-ray-based structure could not provide an answer to the question of whether the double lipid bilayer configuration of AQP0 could functionally maintain a water impermeable state of the channel. To address these questions we conducted molecular dynamics simulations to compare the time-dependent behavior of the AQP0 and AQP1 channels within lipid bilayers. The simulations demonstrate that AQP0, in single or double lipid bilayers, is not closed to water transport and that thermal motions of critical side-chains are sufficient to facilitate the movement of water past any of its constriction regions. These motional requirements do however lead to significant free energy barriers and help explain physiological observations that found water permeability in AQP0 to be substantially lower than in the AQP1 pore.     

Pathogenic mutations in the hydrophobic core of the human prion protein can promote structural instability and misfolding

Transmissible spongiform encephalopathies, or prion diseases, are caused by misfolding and aggregation of the prion protein PrP. These diseases can be hereditary in humans and four of the many disease-associated missense mutants of PrP are in the hydrophobic core: V180I, F198S, V203I and V210I. The T183A mutation is related to the hydrophobic core mutants as it is close to the hydrophobic core and known to cause instability. We used extensive molecular dynamics simulations of these five PrP mutants to compare their dynamics and conformations to those of the wild type PrP. The simulations highlight the changes that occur upon introduction of mutations and help to rationalize experimental findings. Changes can occur around the mutation site, but they can also be propagated over long distances. In particular, the F198S and T183A mutations lead to increased flexibility in parts of the structure that are normally stable, and the short β-sheet moves away from the rest of the protein. Mutations V180I, V210I and, to a lesser extent, V203I cause changes similar to those observed upon lowering the pH, which has been linked to misfolding. Early misfolding is observed in one V180I simulation. Overall, mutations in the hydrophobic core have a significant effect on the dynamics and stability of PrP, including the propensity to misfold, which helps to explain their role in the development of familial prion diseases.