Biomedical subjects

C M Soares

Publications and source records attributed to C M Soares.

At least 19 recordsLinked to original sources

PubMed2001

[NiFe] hydrogenase from Desulfovibrio desulfuricans ATCC 27774: gene sequencing, three-dimensional structure determination and refinement at 1.8 A and modelling studies of its interaction with the tetrahaem cytochrome c3.

The primary and three-dimensional structures of a [NiFe] hydrogenase isolated from D. desulfitricans ATCC 27774 were determined, by nucleotide analysis and single-crystal X-ray crystallography. The three-dimensional structural model was refined to R=0.167 and Rfree=0.223 using data to 1.8 A resolution. Two unique structural features are observed: the [4Fe-4S] cluster nearest the [NiFe] centre has been modified [4Fe-3S-3O] by loss of one sulfur atom and inclusion of three oxygen atoms; a three-fold disorder was observed for Cys536 which binds to the nickel atom in the [NiFe] centre. Also, the bridging sulfur atom that caps the active site was found to have partial occupancy, thus corresponding to a partly activated enzyme. These structural features may have biological relevance. In particular, the two less-populated rotamers of Cys536 may be involved in the activation process of the enzyme, as well as in the catalytic cycle. Molecular modelling studies were carried out on the interaction between this [NiFe] hydrogenase and its physiological partner, the tetrahaem cytochrome c3 from the same organism. The lowest energy docking solutions were found to correspond to an interaction between the haem IV region in tetrahaem cytochrome c3 with the distal [4Fe-4S] cluster in [NiFe] hydrogenase. This interaction should correspond to efficient electron transfer and be physiologically relevant, given the proximity of the two redox centres and the fact that electron transfer decay coupling calculations show high coupling values and a short electron transfer pathway. On the other hand, other docking solutions have been found that, despite showing low electron transfer efficiency, may give clues on possible proton transfer mechanisms between the two molecules.

Amino Acid Sequence↗

PubMed2001

Gene cluster of Rhodothermus marinus high-potential iron-sulfur Protein: oxygen oxidoreductase, a caa(3)-type oxidase belonging to the superfamily of heme-copper oxidases.

The respiratory chain of the thermohalophilic bacterium Rhodothermus marinus contains an oxygen reductase, which uses HiPIP (high potential iron-sulfur protein) as an electron donor. The structural genes encoding the four subunits of this HiPIP:oxygen oxidoreductase were cloned and sequenced. The genes for subunits II, I, III, and IV (named rcoxA to rcoxD) are found in this order and seemed to be organized in an operon of at least five genes with a terminator structure a few nucleotides downstream of rcoxD. Examination of the amino acid sequence of the Rcox subunits shows that the subunits of the R. marinus enzyme have homology to the corresponding subunits of oxidases belonging to the superfamily of heme-copper oxidases. RcoxB has the conserved histidines involved in binding the binuclear center and the low-spin heme. All of the residues proposed to be involved in proton transfer channels are conserved, with the exception of the key glutamate residue of the D-channel (E(278), Paracoccus denitrificans numbering). Analysis of the homology-derived structural model of subunit I shows that the phenol group of a tyrosine (Y) residue and the hydroxyl group of the following serine (S) may functionally substitute the glutamate carboxyl in proton transfer. RcoxA has an additional sequence for heme C binding, after the Cu(A) domain, that is characteristic of caa(3) oxidases belonging to the superfamily. Homology modeling of the structure of this cytochrome domain of subunit II shows no marked electrostatic character, especially around the heme edge region, suggesting that the interaction with a redox partner is not of an electrostatic nature. This observation is analyzed in relation to the electron donor for this caa(3) oxidase, the HiPIP. In conclusion, it is shown that an oxidase, which uses an iron-sulfur protein as an electron donor, is structurally related to the caa(3) class of heme-copper cytochrome c oxidases. The data are discussed in the framework of the evolution of oxidases within the superfamily of heme-copper oxidases.

Amino Acid Sequence↗

PubMed2000

Evolution of the thyroid hormone-binding protein, transthyretin.

Transthyretin (TTR) belongs to a group of proteins, which includes thyroxine-binding globulin and albumin, that bind to and transport thyroid hormones in the blood. TTR is also indirectly implicated in the carriage of vitamin A through the mediation of retinol-binding protein (RBP). It was first identified in 1942 in human serum and cerebrospinal fluid and was formerly called prealbumin for its ability to migrate faster than serum albumin on electrophoresis of whole plasma. It is a single polypeptide chain of 127 amino acids (14,000 Da) and is present in the plasma as a tetramer of noncovalently bound monomers. The major sites of synthesis of TTR in eutherian mammals, marsupials, and birds are the liver and choroid plexus but in reptiles it is synthesised only in the choroid plexus. The observation that TTR is strongly expressed in the choroid plexus but not in the liver of the stumpy-tailed lizard and the strong conservation of expression in the choroid plexus from reptiles to mammals have been taken as evidence to suggest that extrahepatic synthesis of TTR evolved first. The identification and cloning of TTR from the liver of an amphibian, Rana catesbeiana, and a teleost fish, Sparus aurata, and its absence from the choroid plexus of both species suggest an alternative model for its evolution. Protein modelling studies are presented that demonstrate differences in the electrostatic characteristics of the molecule in human, rat, chicken, and fish, which may explain why, in contrast to TTR from human and rat, TTR from fish and birds preferentially binds triiodo-l-thyronine.

Amino Acid Sequence↗

PubMed2000

Ionic strength dependence of the non-physiological electron transfer between flavodoxin and cytochrome c553 from D vulgaris.

A hypothetical model for the non-physiological electron transfer complex between cytochrome c553 (c553) and the flavodoxin (fld) from the sulphate-reducing bacteria Desulfovibrio vulgaris has been recently published [1] based on rigid-body docking and refined by molecular dynamics. In this study, the functional validity of this model is tested by looking at the role of electrostatics in the non-physiological interprotein electron transfer between the two proteins at different ionic strengths. The results are compared with the electron transfer between fld and cytochrome c from horse heart (hhc). Second-order rate constants (k2) were measured for both non-physiological systems at different ionic strengths: a complex, bell-shaped behaviour is observed for the k2 of the c553/fld redox pair with an optimum rate at I=58 mmol l(-1), whereas under the same conditions the k2 for hhc/fld decreased monotonically with increasing ionic strength. Results from the electron transfer kinetics are rationalised in terms of reorganisational effects of an ensemble of conformations of the electron transfer competent c553/fld complexes, consistent with the published model.

Cytochrome c Group↗

PubMed2000

Molecular identification of Paracoccidioides brasiliensis by PCR amplification of ribosomal DNA.

We have amplified and sequenced the 5.8S and 28S ribosomal DNA genes and intergenic regions of Paracoccidioides brasiliensis, strain Pb01. Using primers specifically designed for both ribosomal DNA regions, we were able to discriminate between P. brasiliensis and other human pathogenic fungi by PCR. The use of this molecular marker could be important for paracoccidiodomycosis diagnosis and ecological and molecular epidemiological studies of P. brasiliensis in Latin America.

Base Sequence↗

PubMed2000

Serotypes and subgroups of rotavirus isolated from children in central Brazil.

Group A rotavirus, obtained from children of Goiânia, Brazil, during 1987-1994, were analyzed for subgroup and G serotype by enzyme-linked immunosorbent assay with monoclonal antibodies. The index of serotyping obtained was 61.4% with the following proportions: G1--19.7%, G2--28.0%, G3--9.8%, G4--1.5%, and G5--2.3%. It was observed that G1 occurred from 1987 to 1989 and from 1993 to 1994, and G2 from 1990 to 1993. About 94% of the samples (85/90) could be subgrouped with the following results: 55.5% for SG II, 7.8% SG I, and 31.1% for SG non-I-non-II. Unusual relationship patterns were also detected among serotypes, subgroups, and profiles of electropherotypes in 57.0% of the samples: 20 of them were G2/SG II/"long" profile. The results suggest that variation in temporal and regional characteristics should be considered in the development of rotavirus vaccine.

Antibodies, Monoclonal↗

PubMed2000

Molecular typing of pathogenic fungi.

In this Round Table, the application of several methods of molecular typing were discussed in reference to four important pathogenic fungi: Coccidioides immitis, Histoplasma capsulatum, Candida albicans and Paracoccidioides brasiliensis. Among the different methods the following were discussed: restriction fragment length polymorphisms (RFLP), single nucleotide polymorphisms, random amplified polymorphic DNA (RAPD), polymerase chain reaction (PCR)-RFLP and microsatellites. By means of these methods, several important biological questions related to speciation, mode of reproduction and population genetics could be approached. The basic information obtained from this approach has implications in the understanding of these pathogenic fungi in relation to their behavior and the development of pathogenic features, such as resistance to antimicrobials and virulence factors used for colonization of mammalian hosts. The knowledge obtained from these studies could also be used for the development of innovative diagnostic methods, as well as for novel therapeutic approaches and production of vaccines.

Genetic Techniques↗

PubMed1999

Improved modeling of side-chains in proteins with rotamer-based methods: a flexible rotamer model.

Side-chain modeling has a widespread application in many current methods for protein tertiary structure determination, prediction, and design. Of the existing side-chain modeling methods, rotamer-based methods are the fastest and most efficient. Classically, a rotamer is conceived as a single, rigid conformation of an amino acid sidechain. Here, we present a flexible rotamer model in which a rotamer is a continuous ensemble of conformations that cluster around the classic rigid rotamer. We have developed a thermodynamically based method for calculating effective energies for the flexible rotamer. These energies have a one-to-one correspondence with the potential energies of the rigid rotamer. Therefore, the flexible rotamer model is completely general and may be used with any rotamer-based method in substitution of the rigid rotamer model. We have compared the performance of the flexible and rigid rotamer models with one side-chain modeling method in particular (the self-consistent mean field theory method) on a set of 20 high quality crystallographic protein structures. For the flexible rotamer model, we obtained average predictions of 85.8% for chi1, 76.5% for chi1+2 and 1.34 A for root-mean-square deviation (RMSD); the corresponding values for core residues were 93.0%, 87.7% and 0.70 A, respectively. These values represent improvements of 7.3% for chi1, 8.1% for chi1+2 and 0.23 A for RMSD over the predictions obtained with the rigid rotamer model under otherwise identical conditions; the corresponding improvements for core residues were 6.9%, 10.5% and 0.43 A, respectively. We found that the predictions obtained with the flexible rotamer model were also significantly better than those obtained for the same set of proteins with another state-of-the-art side-chain placement method in the literature, especially for core residues. The flexible rotamer model represents a considerable improvement over the classic rigid rotamer model. It can, therefore, be used with considerable advantage in all rotamer-based methods commonly applied to protein tertiary structure determination, prediction, and design and also in predictions of free energies in mutational studies.

Amino Acids↗

PubMed1999

Crystal structure of cardosin A, a glycosylated and Arg-Gly-Asp-containing aspartic proteinase from the flowers of Cynara cardunculus L.

Aspartic proteinases (AP) have been widely studied within the living world, but so far no plant AP have been structurally characterized. The refined cardosin A crystallographic structure includes two molecules, built up by two glycosylated peptide chains (31 and 15 kDa each). The fold of cardosin A is typical within the AP family. The glycosyl content is described by 19 sugar rings attached to Asn-67 and Asn-257. They are localized on the molecular surface away from the conserved active site and show a new glycan of the plant complex type. A hydrogen bond between Gln-126 and Manbeta4 renders the monosaccharide oxygen O-2 sterically inaccessible to accept a xylosyl residue, therefore explaining the new type of the identified plant glycan. The Arg-Gly-Asp sequence, which has been shown to be involved in recognition of a putative cardosin A receptor, was found in a loop between two beta-strands on the molecular surface opposite the active site cleft. Based on the crystal structure, a possible mechanism whereby cardosin A might be orientated at the cell surface of the style to interact with its putative receptor from pollen is proposed. The biological implications of these findings are also discussed.

Amino Acid Sequence↗

PubMed1999

The caa3 terminal oxidase of the thermohalophilic bacterium Rhodothermus marinus: a HiPIP:oxygen oxidoreductase lacking the key glutamate of the D-channel.

The respiratory chain of the thermohalophilic bacterium Rhodothermus marinus contains a novel complex III and a high potential iron-sulfur protein (HiPIP) as the main electron shuttle (Pereira et al., Biochemistry 38 (1999) 1268-1275 and 1276-1283). In this paper, one of the terminal oxidases expressed in this bacterium is extensively characterised. It is a caa3-type oxidase, isolated with four subunits (apparent molecular masses of 42, 19 and 15 kDa and a C-haem containing subunit of 35 kDa), which has haems of the A(s) type. This oxidase is capable of using TMPD and horse heart cytochrome c as substrates, but has a higher turnover with HiPIP, being the first example of a HiPIP:oxygen oxidoreductase. The oxidase has unusually low reduction potentials of 260 (haem C), 255 (haem A) and 180 mV (haem A3). Subunit I of R. marinus caa3 oxidase has an overall significant homology with the subunits I of the COX type oxidases, namely the metal binding sites and most residues considered to be functionally important for proton uptake and pumping (K- and D-channels). However, a major difference is present: the putative essential glutamate (E278 in Paraccocus denitrificans) of the D-channel is missing in the R. marinus oxidase. Homology modelling of the R. marinus oxidase shows that the phenol group of a tyrosine residue may occupy a similar spatial position as the glutamate carboxyl, in relation to the binuclear centre. Moreover, sequence comparisons reveal that several enzymes lacking that glutamate have a conserved substitution pattern in helix VI: -YSHPXV- instead of -XGHPEV-. These observations are discussed in terms of the mechanisms for proton uptake and it is suggested that, in these enzymes, tyrosine may play the role of the glutamate in the proton channel.

Amino Acid Sequence↗

PubMed1999

An iterative structure-assisted approach to sequence alignment and comparative modeling.

Correct alignment of the sequence of a target protein with those of homologues of known three-dimensional structure is a key step in comparative modeling. Usually an iterative approach that takes account of the local and overall structural features is required. We describe such an approach that exploits databases of structural alignments of homologous proteins (HOMSTRAD, http:/(/)www-cryst.bioc.cam.ac.uk/ approximately homstrad) and protein superfamilies (CAMPASS, http:/(/)www-cryst.bioc.cam.ac.uk/ approximately campass), in which structure-based alignments are analyzed and formatted with the program JOY (http:/(/)www-cryst.bioc.cam.ac.uk/ approximately joy) to reveal conserved local structural features. The databases facilitate the recognition of a family or superfamily, they assist in the selection of useful parent structures, they are helpful in alignment of the target sequences with the parent set, and are useful for deriving relationships that can be used in validating models. In the iterative approach, a model is constructed on the basis of the proposed sequence alignment and this is then reexpressed in the JOY format and realigned with the parent set. This is repeated until the model and sequence alignment is optimized. We examine the case for comparison and use of multiple structures of family members, rather than a single parent structure. We use the targets attempted by our group in CASP3 to assess the value of such procedures.

Algorithms↗

PubMed1999

Improvement of side-chain modeling in proteins with the self-consistent mean field theory method based on an analysis of the factors influencing prediction.

With the objective of improving side-chain conformation prediction, we have analyzed the influence of various factors on prediction by the Self-Consistent Mean Field Theory method, applied to a set of high resolution x-ray protein structure models. These factors may be classed as variations in the mean field optimization protocol, variations in the potential energy function, and variations in rotamer library completeness. We have developed an optimization protocol that consistently reached lower mean field conformational free energies than two other protocols. This protocol led to an important improvement in prediction. We observed a major improvement in prediction with two more detailed van der Waals parameter sets, which we found to be due mainly to the introduction of scaling of 1-4 interactions. In a comparison of two knowledge-based rotamer libraries of considerably different size, we observed an unexpected decrease in prediction with an increase in library completeness. However, when we introduced a torsion potential term in the potential energy function, we found an important increase in average prediction and in the prediction of almost all residue types with a more complete rotamer set. The two knowledge-based rotamer libraries now became equivalent in terms of average prediction. The results we obtained in an analysis of the effect of the introduction of an additional electrostatic term in the potential energy function were largely inconclusive. However, we found a small increase in average prediction for an electrostatic potential term with a fixed dielectric constant of 15. The combined effect of all the factors we analyzed in this study resulted in average prediction accuracies of 79.9% for X1, 68.1% for X1 + 2, and 1.590 A for global rms deviation (RMSD); the corresponding values for core residues were 88.2%, 78.6%, and 1.171 A. These values represent improvements in average prediction of 6.5% for X1, 9.1% for X1 + 2, and 0.163 A for global RMSD over the original conditions; the corresponding improvements in the core were 5.9%, 9.0%, and 0.180 A, respectively.

Amino Acids↗

PubMed1999

Comparative redox and pKa calculations on cytochrome c3 from several Desulfovibrio species using continuum electrostatic methods.

A comparative study of the pH-dependent redox mechanisms of several members of the cytochrome c3 family has been carried out. In a previous work, the molecular determinants of this dependency (the so-called redox-Bohr effect) were investigated for one species using continuum electrostatic methods to find groups with a titrating range and strength of interaction compatible with a mediating role in the redox-Bohr effect. Here we clarify these aspects in the light of new and improved pKa calculations, our findings supporting the hypothesis of propionate D from heme I being the main effector in the pH-dependent modulation of the cytochrome c3 redox potentials in all the c3 molecules studied here. However, the weaker (but significant) role of other titrating groups cannot be excluded, their importance and identity changing with the particular molecule under study. We also calculate the relative redox potentials of the four heme centers among the selected members of the c3 family, using a continuum electrostatic method that takes into account both solvation and interaction effects. Comparison of the calculated values with available data for the microscopic redox potentials was undertaken, the quality of the agreement being dependent upon the choice of the dielectric constant for the protein interior. We find that high dielectric constants give best correlations, while low values result in better magnitudes for the calculated potentials. The possibility that the crystallographic calcium ion in c3 from Desulfovibrio gigas may be present in the solution structure was tested, and found to be likely.

Calcium↗

PubMed1999

Effects of protein-protein interactions on electron transfer: docking and electron transfer calculations for complexes between flavodoxin and c-type cytochromes.

Theoretical studies of protein-protein association and electron transfer were performed on the binary systems formed by Desulfovibrio vulgaris Hildenborough (D. v. H.) flavodoxin and D. v. H. cytochrome c553 and by flavodoxin and horse heart cytochrome c. Initial structures for the complexes were obtained by rigid-body docking and were refined by MD to allow for molecular flexibility. The structures thus obtained were analysed in terms of their relative stability through the calculation of excess energies. Electrostatic, van der Waals and solvation energy terms showed all to have significant contributions to the stability of complexes. In the best association solutions found for both cytochromes, these bind to different zones of flavodoxin. The binding site of flavodoxin observed for cytochrome c is in accordance with earlier works [27]. The various association modes found were characterised in terms of electron transfer using the Pathways model. For complexes between flavodoxin and horse heart cytochrome c, some correlation was observed between electron tunnelling coupling factors and conformation energy; the best conformation found for electron transfer corresponded also to the best one in terms of energy. For complexes between flavodoxin and cytochrome c553 this was not the case and a lower correlation was observed between electron tunnelling coupling factors and excess energies. These results are in accordance with the differences in the experimental dependence of electron transfer rates with ionic strength observed between these two cases.

Cytochrome c Group↗

PubMed1999

Nine-haem cytochrome c from Desulfovibrio desulfuricans ATCC 27774:primary sequence determination, crystallographic refinement at 1.8 and modelling studies of its interaction with the tetrahaem cytochrome c3.

A monomeric nine-haem cytochrome c (9Hcc) with 292 amino acid residues was isolated from cells of the sulfate- and nitrate-reducing bacterium Desulfovibrio desulfuricans ATCC 27774 grown under both nitrate- and sulfate-respiring conditions. The nucleotide sequence encoding the 292 residues was determined, allowing the correction of about 10% of the previous primary structure, determined from 1.8 A electron density maps. The refinement at 1.8 A resolution of the structural model was completed, giving an R-value of 16.5%. The nine haem groups are arranged into two tetrahaem clusters, located at both ends of the molecule, with Fe-Fe distances and local protein fold very similar to tetrahaem cytochromes c3, and the extra haem is located asymmetrically between the two regions. The new primary sequence determination confirmed the 39% sequence homology found between this cytochrome and the C-terminal region (residues 229-514) of the high-molecular-weight cytochrome c (Hmc) from D. vulgaris Hildenborough, providing strong evidence of structural similarity between 9Hcc and the C-terminal region of Hmc. The interaction between 9Hcc and the tetrahaem cytochrome c3 from the same organism was studied by modelling methods, and the results suggest that a specific interaction is possible between haem 4 of tetrahaem cytochrome c3 and haem 1 or haem 2 of 9Hcc, in agreement with previous kinetic experiments which showed the catalytic effect of the tetrahaem cytochrome c3 upon the reduction of 9Hcc by the [NiFe] hydrogenase from D. desulfuricans ATCC 27774. These studies suggest a role for 9Hcc as part of the assembly of redox proteins involved in recycling the molecular hydrogen released by the cell as a result of substrate oxidation.

Amino Acid Sequence↗

PubMed1999

Simulation of electron-proton coupling with a Monte Carlo method: application to cytochrome c3 using continuum electrostatics.

A new method is presented for simulating the simultaneous binding equilibrium of electrons and protons on protein molecules, which makes it possible to study the full equilibrium thermodynamics of redox and protonation processes, including electron-proton coupling. The simulations using this method reflect directly the pH and electrostatic potential of the environment, thus providing a much closer and realistic connection with experimental parameters than do usual methods. By ignoring the full binding equilibrium, calculations usually overlook the twofold effect that binding fluctuations have on the behavior of redox proteins: first, they affect the energy of the system by creating partially occupied sites; second, they affect its entropy by introducing an additional empty/occupied site disorder (here named occupational entropy). The proposed method is applied to cytochrome c3 of Desulfovibrio vulgaris Hildenborough to study its redox properties and electron-proton coupling (redox-Bohr effect), using a continuum electrostatic method based on the linear Poisson-Boltzmann equation. Unlike previous studies using other methods, the full reduction order of the four hemes at physiological pH is successfully predicted. The sites more strongly involved in the redox-Bohr effect are identified by analysis of their titration curves/surfaces and the shifts of their midpoint redox potentials and pKa values. Site-site couplings are analyzed using statistical correlations, a method much more realistic than the usual analysis based on direct interactions. The site found to be more strongly involved in the redox-Bohr effect is propionate D of heme I, in agreement with previous studies; other likely candidates are His67, the N-terminus, and propionate D of heme IV. Even though the present study is limited to equilibrium conditions, the possible role of binding fluctuations in the concerted transfer of protons and electrons under nonequilibrium conditions is also discussed. The occupational entropy contributions to midpoint redox potentials and pKa values are computed and shown to be significant.

Binding Sites↗

PubMed1999

Differential expression of an hsp70 gene during transition from the mycelial to the infective yeast form of the human pathogenic fungus Paracoccidioides brasiliensis.

We have isolated and characterized cDNA and genomic clones that encode a 70 kDa heat shock protein (Hsp70) from the dimorphic human pathogenic fungus Paracoccidioides brasiliensis. The gene encodes a 649-amino-acid protein showing high identity with other members of the hsp70 gene family. The hsp70 gene is induced during both heat shock of yeast cells at 42 degrees C and the mycelial to yeast transition. A differential expression of this gene can be observed between mycelial and yeast forms, with a much higher level of expression in the yeast. We found two introns of 178 and 72 nucleotides in the P. brasiliensis hsp70 gene. Splicing of these introns is regulated during the heat shock process and possibly during infection. In order to analyse the differential accumulation of unspliced mRNA following cellular differentiation and/or heat shock, reverse transcriptase-polymerase chain reaction (RT-PCR) experiments were carried out. The temperature-induced mycelial to yeast transition results in the transient accumulation of unspliced hsp70 mRNA transcripts. Yeast cells, after adaptation at 36 degrees C, seem to be more proficient at splicing, at least with respect to hsp70 mRNA because, during a severe heat shock (42 degrees C), the unspliced form of this mRNA does not accumulate. The mycelial to yeast differentiation will have the adaptational effect of increasing the resistance of the organism to environmental stress, which may be necessary for parasite survival in the mammalian host.

Amino Acid Sequence↗

PubMed1999

Local water bridges and protein conformational stability.

Recent studies have pointed out the important role of local water structures in protein conformational stability. Here, we present an accurate and computationally effective way to estimate the free energy contribution of the simplest water structure motif--the water bridge. Based on the combination of empirical parameters for accessible protein surface area and the explicit consideration of all possible water bridges with the protein, we introduce an improved protein solvation model. We find that accounting for water bridge formation in our model is essential to understand the conformational behavior of polypeptides in water. The model formulation, in fact, does not depend on the polypeptide nature of the solute and is therefore applicable to other flexible biomolecules (i.e., DNAs, RNAs, polysaccharides, etc.).

Hydrogen Bonding↗