IGHD II: A Novel GH-1 Gene Mutation (GH-L76P) Severely Affects GH Folding, Stability, and Secretion

CONTEXT The autosomal dominant form of GH deficiency (IGHD II) is characterized by markedly reduced GH secretion combined with low concentrations of IGF-1 leading to short stature. OBJECTIVE Structure-function analysis of a missense mutation in the GH-1 gene converting codon 76 from leucine (L) to proline (P) yielding a mutant GH-L76P peptide. DESIGN, SETTINGS, AND PATIENTS Heterozygosity for GH-L76P/wt-GH was identified in a nonconsanguineous Spanish family. The index patients, two siblings, a boy and a girl, were referred for assessment of their short stature (-3.2 and -3.8 SD). Their grandmother, father, and aunt were also carrying the same mutation and showed severe short stature; therefore, IGHD II was diagnosed. INTERVENTIONS AND RESULTS AtT-20 cells coexpressing both wt-GH and GH-L76P showed a reduced GH secretion (P < .001) after forskolin stimulation when compared with the cells expressing only wt-GH. In silico mutagenesis and molecular dynamics simulations presented alterations of correct folding and mutant stability compared with wt-GH. Therefore, further structural analysis of the GH-L76P mutant was performed using expressed and purified proteins in Escherichia coli by thermofluor assay and fast degradation proteolysis assay. Both assays revealed that the GH-L76P mutant is unstable and misfolded compared to wt-GH confirming the bioinformatic model prediction. CONCLUSIONS This is the first report of a family suffering from short stature caused by IGHD II, which severely affects intracellular GH folding and stability as well as secretion, highlighting the necessity of functional analysis of any GH variant for defining new mechanisms as a cause for IGHD II.

Targeting TRAF6 for Cancer Therapeutical Development

Tumor necrosis factor (TNF)-Receptor Associated Factors (TRAFs) are a family of signal transducer proteins. TRAF6 is a unique member of this family in that it is involved in not only the TNF superfamily, but the toll-like receptor (TLR)/IL-1R (TIR) superfamily. The formation of the complex consisting of Receptor Activator of Nuclear Factor κ B (RANK), with its ligand (RANKL) results in the recruitment of TRAF6, which activates NF-κB, JNK and MAP kinase pathways. TRAF6 is critical in signaling with leading to release of various growth factors in bone, and promotes osteoclastogenesis. TRAF6 has also been implicated as an oncogene in lung cancer and as a target in multiple myeloma. In the hopes of developing small molecule inhibitors of the TRAF6-RANK interaction, multiple steps were carried out. Computational prediction of hot spot residues on the protein-protein interaction of TRAF6 and RANK were examined. Three methods were used: Robetta, KFC2, and HotPoint, each of which uses a different methodology to determine if a residue is a hot spot. These hot spot predictions were considered the basis for resolving the binding site for in silico high-throughput screening using GOLD and the MyriaScreen database of drug/lead-like compounds. Computationally intensive molecular dynamics simulations highlighted the binding mechanism and TRAF6 structural changes upon hit binding. Compounds identified as hits were verified using a GST-pull down assay, comparing inhibition to a RANK decoy peptide. Since many drugs fail due to lack of efficacy and toxicity, predictive models for the evaluation of the LD50 and bioavailability of our TRAF6 hits, and these models can be used towards other drugs and small molecule therapeutics as well. Datasets of compounds and their corresponding bioavailability and LD50 values were curated based, and QSAR models were built using molecular descriptors of these compounds using the k-nearest neighbor (k-NN) method, and quality of these models were cross-validated.

Li & Pb MD Sinulations

Pb17Li is today a reference breeder material in diverse fusion R&D programs worldwide. Extracting dynamic and structural properties of liquid LiPb mixtures via molecular dynamics simulations, represent a crucial step for multiscale modeling efforts in order to understand the suitability of this compound for future Nuclear Fusion technologies. At present a Li-Pb cross potential is not available in the literature. Here we present our first results on the validation of two semi-empirical potentials for Li and Pb in liquid phase. Our results represent the establishment of a solid base as a previous but crucial step to implement a LiPb cross potential. Structural and thermodynamical analyses confirm that the implemented potentials for Li and Pb are realistic to simulate both elements in the liquid phase.

Assessment of interatomic potentials for atomistic analysis of static and dynamic properties of screw dislocations in W

Screw dislocations in bcc metals display non-planar cores at zero temperature which result in high lattice friction and thermally-activated strain rate behavior. In bcc W, electronic structure molecular statics calculations reveal a compact, non-degenerate core with an associated Peierls stress between 1.7 and 2.8 GPa. However, a full picture of the dynamic behavior of dislocations can only be gained by using more efficient atomistic simulations based on semiempirical interatomic potentials. In this paper we assess the suitability of five different potentials in terms of static properties relevant to screw dislocations in pure W. Moreover, we perform molecular dynamics simulations of stress-assisted glide using all five potentials to study the dynamic behavior of screw dislocations under shear stress. Dislocations are seen to display thermally-activated motion in most of the applied stress range, with a gradual transition to a viscous damping regime at high stresses. We find that one potential predicts a core transformation from compact to dissociated at finite temperature that affects the energetics of kink-pair production and impacts the mechanism of motion. We conclude that a modified embedded-atom potential achieves the best compromise in terms of static and dynamic screw dislocation properties, although at an expense of about ten-fold compared to central potentials.

Representing microstates of static granular matter in a stress phase space

A stress phase space is proposed to compare the static packings of a granular system (microstates) that are compatible to a macrostate described by external stresses. The equivalent stress of each particle of a static packing can be obtained from the mechanical interaction forces, and the associated volume is given by the respective Voronoi cell. Therefore, particles can be located at different stress levels and grouped into categories or configurations, which are defined in base of the geometrical features of the local arrangement (in particular, of the number of forces that keep them force-balanced). They can be represented as points in a stress phase space. The nature of this space is analyzed in detail. The integration limits of the stress variables that avoid or limit tensile states and the capability of each configuration to represent specific stress states establish its main features. Furthermore, if some stress variables are used, instead of the usual components of the Cauchy stress tensor, then some symmetries can be found. Results obtained from molecular dynamics simulations are used to check this nature. Finally, some statistical ensembles are written in terms of the coordinates of this phase space. These require some assumptions that are made in base on continuum mechanics principles.

The role of enzyme distortion in the single displacement mechanism of family 19 chitinases

By using molecular dynamics simulations, we have examined the binding of a hexaNAG substrate and two potential hydrolysis intermediates (an oxazoline ion and an oxocarbenium ion) to a family 19 barley chitinase. We find the hexaNAG substrate binds with all sugars in a chair conformation, unlike the family 18 chitinase which causes substrate distortion. Glu 67 is in a position to protonate the anomeric oxygen linking sugar residues D and E whereas Asn 199 serves to hydrogen bond with the C2′ N-acetyl group of sugar D, thus preventing the formation of an oxazoline ion intermediate. In addition, Glu 89 is part of a flexible loop region allowing a conformational change to occur within the active site to bring the oxocarbenium ion intermediate and Glu 89 closer by 4–5 Å. A hydrolysis product with inversion of the anomeric configuration occurs because of nucleophilic attack by a water molecule that is coordinated by Glu 89 and Ser 120. Issues important for the design of inhibitors specific to family 19 chitinases over family 18 chitinases also are discussed.

Ionophores and receptors using cation-π interactions: Collarenes

Cation-π interactions are important forces in molecular recognition by biological receptors, enzyme catalysis, and crystal engineering. We have harnessed these interactions in designing molecular systems with circular arrangement of benzene units that are capable of acting as ionophores and models for biological receptors. [n]Collarenes are promising candidates with high selectivity for a specific cation, depending on n, because of their structural rigidity and well-defined cavity size. The interaction energies of [n]collarenes with cations have been evaluated by using ab initio calculations. The selectivity of these [n]collarenes in aqueous solution was revealed by using statistical perturbation theory in conjunction with Monte Carlo and molecular dynamics simulations. It has been observed that in [n]collarenes the ratio of the interaction energies of a cation with it and the cation with the basic building unit (benzene) can be correlated to its ion selectivity. We find that collarenes are excellent and efficient ionophores that bind cations through cation-π interactions. [6]Collarene is found to be a selective host for Li+ and Mg2+, [8]collarene for K+ and Sr2+, and [10]collarene for Cs+ and Ba2+. This finding indicates that [10]collarene and [8]collarene could be used for effective separation of highly radioactive isotopes, 137Cs and 90Sr, which are major constituents of nuclear wastes. More interestingly, collarenes of larger cavity size can be useful in capturing organic cations. [12]Collarene exhibits a pronounced affinity for tetramethylammonium cation and acetylcholine, which implies that it could serve as a model for acetylcholinestrase. Thus, collarenes can prove to be novel and effective ionophores/model-receptors capable of heralding a new direction in molecular recognition and host-guest chemistry.

A molecular level picture of the stabilization of A-DNA in mixed ethanol–water solutions

Advances in computer power, methodology, and empirical force fields now allow routine “stable” nanosecond-length molecular dynamics simulations of DNA in water. The accurate representation of environmental influences on structure remains a major, unresolved issue. In contrast to simulations of A-DNA in water (where an A-DNA to B-DNA transition is observed) and in pure ethanol (where disruption of the structure is observed), A-DNA in ≈85% ethanol solution remains in a canonical A-DNA geometry as expected. The stabilization of A-DNA by ethanol is likely due to disruption of the spine of hydration in the minor groove and the presence of ion-mediated interhelical bonds and extensive hydration across the major groove.

Mapping the early steps in the pH-induced conformational conversion of the prion protein

Under certain conditions, the prion protein (PrP) undergoes a conformational change from the normal cellular isoform, PrPC, to PrPSc, an infectious isoform capable of causing neurodegenerative diseases in many mammals. Conversion can be triggered by low pH, and in vivo this appears to take place in an endocytic pathway and/or caveolae-like domains. It has thus far been impossible to characterize the conformational change at high resolution by experimental methods. Therefore, to investigate the effect of acidic pH on PrP conformation, we have performed 10-ns molecular dynamics simulations of PrPC in water at neutral and low pH. The core of the protein is well maintained at neutral pH. At low pH, however, the protein is more dynamic, and the sheet-like structure increases both by lengthening of the native β-sheet and by addition of a portion of the N terminus to widen the sheet by another two strands. The side chain of Met-129, a polymorphic codon in humans associated with variant Creutzfeldt–Jakob disease, pulls the N terminus into the sheet. Neutralization of Asp-178 at low pH removes interactions that inhibit conversion, which is consistent with the Asp-178–Asn mutation causing human prion diseases.

Protein folding from a highly disordered denatured state: The folding pathway of chymotrypsin inhibitor 2 at atomic resolution

Previous experimental and theoretical studies have produced high-resolution descriptions of the native and folding transition states of chymotrypsin inhibitor 2 (CI2). In similar fashion, here we use a combination of NMR experiments and molecular dynamics simulations to examine the conformations populated by CI2 in the denatured state. The denatured state is highly unfolded, but there is some residual native helical structure along with hydrophobic clustering in the center of the chain. The lack of persistent nonnative structure in the denatured state reduces barriers that must be overcome, leading to fast folding through a nucleation–condensation mechanism. With the characterization of the denatured state, we have now completed our description of the folding/unfolding pathway of CI2 at atomic resolution.

Comparison of the early stages of forced unfolding for fibronectin type III modules

The structural changes accompanying stretch-induced early unfolding events were investigated for the four type III fibronectin (FN-III) modules, FN-III7, FN-III8, FN-III9, and FN-III10 by using steered molecular dynamics. Simulations revealed that two main energy barriers, I and II, have to be overcome to initiate unraveling of FN-III's tertiary structure. In crossing the first barrier, the two opposing β-sheets of FN-III are rotated against each other such that the β-strands of both β-sheets align parallel to the force vector (aligned state). All further events in the unfolding pathway proceed from this intermediate state. A second energy barrier has to be overcome to break the first major cluster of hydrogen bonds between adjacent β-strands. Simulations revealed that the height of barrier I varied significantly among the four modules studied, being largest for FN-III7 and lowest for FN-III10, whereas the height of barrier II showed little variation. Key residues affecting the mechanical stability of FN-III modules were identified. These results suggest that FN-III modules can be prestretched into an intermediate state with only minor changes to their tertiary structures. FN-III10, for example, extends 12 Å from the native “twisted” to the intermediate aligned state, and an additional 10 Å from the aligned state to further unfolding where the first β-strand is peeled away. The implications of the existence of intermediate states regarding the elasticity of fibrillar fibers and the stretch-induced exposure of cryptic sites are discussed.

The arginine finger of RasGAP helps Gln-61 align the nucleophilic water in GAP-stimulated hydrolysis of GTP

The Ras family of GTPases is a collection of molecular switches that link receptors on the plasma membrane to signaling pathways that regulate cell proliferation and differentiation. The accessory GTPase-activating proteins (GAPs) negatively regulate the cell signaling by increasing the slow intrinsic GTP to GDP hydrolysis rate of Ras. Mutants of Ras are found in 25–30% of human tumors. The most dramatic property of these mutants is their insensitivity to the negative regulatory action of GAPs. All known oncogenic mutants of Ras map to a small subset of amino acids. Gln-61 is particularly important because virtually all mutations of this residue eliminate sensitivity to GAPs. Despite its obvious importance for carcinogenesis, the role of Gln-61 in the GAP-stimulated GTPase activity of Ras has remained a mystery. Our molecular dynamics simulations of the p21ras–p120GAP–GTP complex suggest that the local structure around the catalytic region can be different from that revealed by the x-ray crystal structure. We find that the carbonyl oxygen on the backbone of the arginine finger supplied in trans by p120GAP (Arg-789) interacts with a water molecule in the active site that is forming a bridge between the NH2 group of the Gln-61 and the γ-phosphate of GTP. Thus, Arg-789 may play a dual role in generating the nucleophile as well as stabilizing the transition state for P—O bond cleavage.

Direct evidence for modified solvent structure within the hydration shell of a hydrophobic amino acid.

Neutron scattering experiments are used to determine scattering profiles for aqueous solutions of hydrophobic and hydrophilic amino acid analogs. Solutions of hydrophobic solutes show a shift in the main diffraction peak to smaller angle as compared with pure water, whereas solutions of hydrophilic solutes do not. The same difference for solutions of hydrophobic and hydrophilic side chains is also predicted by molecular dynamics simulations. The neutron scattering curves of aqueous solutions of hydrophobic amino acids at room temperature are qualitatively similar to differences between the liquid molecular structure functions measured for ambient and supercooled water. The nonpolar solute-induced expansion of water structure reported here is also complementary to recent neutron experiments where compression of aqueous solvent structure has been observed at high salt concentration.

Mechanical annealing of metallic electrodes at the atomic scale

The process of creating an atomically defined and robust metallic tip is described and quantified using measurements of contact conductance between gold electrodes and numerical simulations. Our experiments show how the same conductance behavior can be obtained for hundreds of cycles of formation and rupture of the nanocontact by limiting the indentation depth between the two electrodes up to a conductance value of approximately 5G0 in the case of gold. This phenomenon is rationalized using molecular dynamics simulations together with density functional theory transport calculations which show how, after repeated indentations (mechanical annealing), the two metallic electrodes are shaped into tips of reproducible structure. These results provide a crucial insight into fundamental aspects relevant to nanotribology or scanning probe microscopies.

Mechanical, electrical, and magnetic properties of Ni nanocontacts

The dynamic deformation upon stretching of Ni nanowires as those formed with mechanically controllable break junctions or with a scanning tunneling microscope is studied both experimentally and theoretically. Molecular dynamics simulations of the breaking process are performed. In addition, and in order to compare with experiments, we also compute the transport properties in the last stages before failure using the first-principles implementation of Landauer's formalism included in our transport package ALACANT.